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The purpose of this research was to investigate the effect of tidal - 
period internal wave vertical motion on oceanic mixed layer dynamics, 
and to discern the effect upon the diurnal evolution of thermoclines, 
as observed during MILE (mixed layer experiment). Vertical advection 
was added to a one-dimensional bulk model of the mixed layer with an 
assumed linear in z and sinusoidal in time dependence. The rate of 
mixed layer deepening was therefore due to the combination of vertical 
mction and entrainment. The first significant result was the finding 
that the interaction between vertical mixing and vertical motion 
depended upon the wave frequency and its phase relation to the diurnal 
heating cycle. Second, linear and non-linear interactions of the wave 
induced vertical motion with the cyclical boundary conditions can 
generate two-dimensional (x-z) structure in the near-surface temperature 
field of an initially horizontally homogeneous ocean under the influence 
of horizontally homogeneous surface boundary conditions. Finally, this 
advective interaction increases the utility of the mixed layer model in 
single station forecasting. 



4 



TABLE OF CONTENTS 



I. INTRODUCTION - 11 

A. BACKGROUND DATA ANALYSIS 12 

B. ADVECTION IN MODELLING 14 

C. SCOPE OF THIS STUDY 15 

II. MIXED LAYER MODEL AND VERTICAL MOTION THEORY 17 

A. MIXED LAYER MODEL 17 

1. Boundary Layer Deepening by Entrainment 20 

2. Boundary Layer Shallowing 21 

B. VERTICAL MOTION THEORY 22 

III. EXPERIMENTS - - 24 

A. PRODUCTS 24 

B. STANDARDS 26 

% 

1 . Standard 1 (SI ) 27 

2. Standard 2 (S2) 30 

C. EXPERIMENTS — 34 

1. Case I — 36 

a. Constant Downwel 1 ing (Case la) 37 

b. Constant Downwel ling (Case lb) 40 

c. Constant Upwelling (Case Ic) 43 

d. Constant Upwelling (Case Id) 45 

2. Case II 48 

a. Pulsed Downwel 1 ing (Case I la) 48 

(1) 180° Out of Phase - 50 



5 



(2) Maximum Downwelling Before Heating 50 

(3) Coupled Heating and Downwelling 51 

b. Pulsed Downwelling (Case I I b ) 53 

(1) 180° Out of Phase With Heating 53 

(2) Maximum Downwelling Before Heating 53 

(3) Coupled Heating and Downwelling — - — 54 

c. Pulsed Upwelling (Case lie) 54 

(1) 180° Out of Phase With Heating 55 

(2) Maximum Upwelling Before Heating 55 

(3) Coupled Heating and Upwelling 56 

d. (Pulsed Upwelling (Case I Id) 56 

(1) 180° Out of Phase With Heating 56 

(2) Maximum Upwelling Before Heating 57 

(3) Coupled Heating And Upwelling 58 

3. Case III 58 

a. Shallow Initial Layer (Case Ilia) — — - 60 

(1) Peak Upward At Hour 7 60 

(2) Peak Upward At Hour 1 61 

(3) Peak Upward Motion At Hour 19 62 

(4) Peak Upward Motion At Hour 13 62 

(5) Two-Dimensional Data 63 

b. Deep Initial Layer (Case Illb) 64 

(1) Peak Upward Motion At Hour 7 65 

(2) Peak Upward Motion At Hour 1 66 

(3) Peak Upward Motion At Hour 19 66 

(4) Peak Upward Motion At Hour 13 — 67 

(5) Two-Dimensional Data 68 



6 



4. Case IV 70 

a. Diurnal Wave (Case IVa) 71 

0) Two-Dimensional Data 71 

b. Semi-diurnal Wave (Case IVb) 73 

0) Six Hour Shift 73 

(2) Twelve Hour Shift 75 

(3) . Two-Dimensional Data 77 

5. Case V (Superimposed Waves ) 78 

IV. CONCLUSIONS AND RECOMMENDATIONS — 82 

FIGURES 11 THROUGH 53 34 

BIBLIOGRAPHY 127 

INITIAL DISTRIBUTION LIST 129 



7 



LIST OF FIGURES 



1 . Hourly parameters, temperature 

contours and profiles for Case SI 25 

2. Initial temperature profile for shallow layer condition — 28 

3. Initial temperature profile for deep layer condition 31 

4. Temperature profile variation from 

t] to t .2 under upwelling conditions 38 

5. Pulsed upwelling form 49 

6. Vertical motion-wave form relationship 59 

7. Diurnal wave relationship to model location 72 

8. Phase relationship of Case IVb (six hour shift) 74 

9. Phase relationship of Case IVb (twelve hour shift) 76 

10. MILE spectral analysis 79 

11. Same as Figure 1 for Case S2 84 

12. As Figure 1, Case la 85 

13. As Figure 1, Case lb 86 

14. As Figure 1, Case Ic 87 

15. As Figure 1, Case Id 88 

16. As Figure 1, Case Ha, peak motion at hour 18 89 

17. As Figure 1, Case Ila, peak motion at hour 24 90 

18. As Figure 1, Case Ila, peak motion at hour 5 91 

19. As Figure 1, Case lib, peak motion at hour 18 92 

20. As Figure 1, Case 1 1 b , peak motion at hour 24 93 

21. As Figure 1, Case 1 1 b , peak motion at hour 5 94 

22. As Figure 1, Case lie, peak motion at hour 24 95 

23. As Figure 1, Case lie, peak motion at hour 5 96 

24. As Figure 1, Case I Id , peak motion at hour 24 97 



8 



25. As Figure 1, Case lid, peak motion at hour 5 98 

26. Case lie and I Id, peak motion at hour 18 99 

27. As Figure 1, Case Ilia, peaks at hours 7 and 19 100 

28. As Figure 1, Case Ilia, peaks at hours 1 and 13 101 

29. As Figure 27, reversed direction of the peak motion 102 

30. As Figure 28, reversed direction of the peak motion 103 

31. Contours for Case Ilia in two dimensional mode 104 

32. Profiles for hours 6 and 12 

in two dimensional mode. Case Ilia 105 

33. Profiles for hours 18 and 24 

in two-dimensional mode. Case Ilia 106 

34. As Figure 1, Case 1 1 1 b , peaks at hours 7 and 19 107 

35. As Figure 1, Case 1 1 1 b , peaks at hours 1 and 13 108 

36. As Figure 34, reversed direction of the peak motion 109 

37. As Figure 35, reversed direction of the peak motion 110 

38. As Figure 31, Case 1 1 1 b 111 

39. As Figure 32, Case 1 1 1 b 112 

40. As Figure 33, Case 1 1 1 b 113 

41. Contours of Case IVa with 4 phases of initial wave 114 

42. Profiles for Case IVa 115 

43. Profiles for Case IVa 116 

44. As Figure 31, Case IVa 117 

45. As Figure 32, Case IVa — 118 

46. As Figure 33, Case IVa 119 

47. As Figure 1, Case IVb, 2 nhase shifts 120 

48. As Figure 31, Case IVb 121 

49. As Figure 32, Case IVb 122 



9 



123 



50. As Figure 33, Case IVb 

51. As Figure 31, Case Va 124 

52. As Figure 32, Case Va 125 

53. As Figure 33, Case Vb 126 



10 



I. INTRODUCTION 



In recent years much study of internal waves has focused on the 
linear and non-linear wave interactions. Models for wave prediction 
have been developed and tested to match the theory. There have been 
attempts to measure the internal wave fields in situ. Some of the 
most useful of these have been the observations of temperature and 
pressure field variations with depth, providing a basis for computa- 
tions of the vertical motions and interactions involved. These studies 
and analyses of data demonstrate the importance of the internal wave 
field in upper ocean and near surface dynamics. 

Roth et al . (1981) presented clear evidence that internal waves 
are not as critical in an analysis of the dynamics of mixing at depth 
because they have little variability. However, near the surface, 
where they interact with other phenomena such as cyclical heat flux, 
wind stress, and other boundary conditions, their effect can be 
considerable when dealing with the budgets of heat and momentum. 
Holloway (1980) pointed out that internal waves are not weak waves 
(in the sense of weak resonant interaction theory) and have a high 
impact on mixed layer dynamics down to and including turbulent scales 
of motion. Pinkel (1981) showed that the near-surface wave field has 
a profound effect on the vertical velocity field in the upper ocean. 

Some mixed layer experiments such as those by Linden (1975) and 
Kantha et al . (1977) have shown that internal wave radiation may play 
an important role in controlling the rate of deepening of a turbulent 
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mixed layer in a stratified fluid. Displacements induced by the turbu- 
lence in the base of the mixed layer served to form internal waves 
which removed kinetic energy of motion from the mixed layer. Since 
this energy is also important in the entrainment process, these induced 
internal waves, of relatively high frequency, serve to retard deepening 
of a mixed layer. Townsend (1968) has derived an expression for the 
rate of generation of the internal waves and their energy by a steady 
disturbance field which is advected relative to a stratified layer at 
a constant speed. The theory was extended by Bell (1978) to cover 
general space-time dependence of the disturbance field. These works 
indicate that internal waves could be a means of retarding the growth 
of mixed layers, through the removal of the kinetic energy in the layers. 

Such studies all indicate that the presence of internal waves should 
be taken into account as a possibly important factor in the mechanical 
energy budget which controls mixing. It is hypothesized that internal 
waves which have periods comparable to the dominant time scales of the 
atmospheric forcing (surface heat flux and wind stress) will be the 
most likely to result in an interaction with the entrainment mixing 
process in terms of a positive feedback (increased mixing) effect. 

A. BACKGROUND DATA ANALYSIS 

During the period August-September 1977, a mixed layer experiment 
(MILE) was conducted at Ocean Station P (55° N, 145° W) in the North 
Pacific. The depth of the mixed layer varied between about five and 
thirty meters during the period of the experiment. Temperature was 
continuously observed at fixed depths for a period of three weeks. A 
spectrum analysis of this data (See Figure 10 in later section) shows 
some important features. 
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First, a major spectral peak was at about 12.79 t .02 hours, with 
secondary peaks at 6,62 + .13 and 3.94 + .06 hours. Lower order 
spectral peaks are at 4.9 hours, and at several frequencies having 
periods under 3.5 hours. The peak at 12.79 hours is clearly asso- 
ciated with the semi-diurnal tide (12.45 hrs.), with the other two 
peaks at 6.62 and 3.94 hours probably being higher harmonics. At all 
three depths of the analysis it is clear that the semi-diurnal internal 
tidal wave, with its harmonics, is the dominant cause of vertical 
motion at the base of the mixed layer. 

Second, the tidal peaks, including the harmonics, are very coherent 
vertically throughout the depths analyzed. They are present in all 
three depths analyzed, below the layer, and in the mixed layer, with 
very little frequency or phase shift. 

Finally, the spectral peaks at other than tidal frequencies are 
not vertically coherent. Also, there is a large variation in non-tidal 
spectral densities for observations above and below the thirty meter 
depth— the top of the seasonal thermocline at this time of year. This 
suggests that something has happened to the higher frequency waves 
near the surface making the interpretation of their effect on mixed 
layer dynamics hard to judge from simple theories. However, at the 
same time it suggests that the effect of tidal internal waves may be 
quite easily discerned. 

These results confirm the work of Pinkel (1981) whose studies 
showed that of the wave fields in the upper ocean, the tide and its 
harmonics have the greatest vertical coherence. 

Studies by Simpson and Paulson (1979) demonstrated the variability 
of the temperature structure in the upper ocean during the POLE experiment. 
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Shook (1980) showed that these variations, were a response to vertical 
and horizontal advective processes. In a study of Pacific Ocean ship 
track temperature data, Butler (1981) found definite cyclical varia- 
tions in the sea surface temperature structure. 

From these studies of actual data, it is obvious that (1) the tidal 
period internal waves and their harmonics are an important source of 
energy in the upper ocean, and (2) surface layer temperature patch- 
iness exists in the ocean, and may in part be generated by interaction 
of mixing and internal wave motion. 

B. ADVECTION IN MODELLING 

To date little has been published concerning the interaction of 
tidal-period internal waves and entrainment by the mixed layer. Ge 
Szoeke (1980) has discussed advection in the mixed layer as a result 
of horizontal variability of the surface wind stress. He showed that 
advection, upwelling and downwelling, was an important parameter in 
mixed layer dynamics, resulting from the boundary forcing of the wind 
stress. Cushman-Roisin (1981), Garwood et al . (1981), and Price (1981) 
also demonstrated that advection can play an important role in the 
dynamics of the mixed layer, but the interaction with internal gravity 
waves is not treated. De Szoeke and Rhines (1981) discussed coastal 
upwelling at length, but the advective motion was conceived as a 
necessary factor to achieve a realistic balance for surface heating, 
transport, and mixing, in order to produce observed coastal temperature 
fields. Price et al . (1978) included vertical advection in the test 
of a mixed layer formulation during storm conditions. They found that 
the phase dependence between wind stress and mixed layer deepening, 
and the mixed layer deepening itself were reasonably simulated. 
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Within this discussion of a limited amount of relevant literature 
on the advective process it has been shown that vertical advection 
must play an important role in storm response, coastal circulation, 
and fronts. No one has examined the role of internal wave motion, on 
the tidal scale, in the open ocean where the geometry and boundary 
conditions are rather one-dimensional. 

In possibly the most relevant previous effort, Spigel (1980) 
looked at the coupling of internal wave motion and entrainment in a 
two-layer lake model, having a mixed layer which is based upon the 
one-dimensional turbulent kinetic energy budget. Time scales of 
internal wave motion, wave decay, and entrainment were widely 
separated. Using a seiche effect in a basin, he demonstrated that 
internal waves did in fact play an important role in the entrainment 
zone. 

With this background it was believed that the time has come for 
an investigation of the effects of internal vertical motion in the 
oceanic mixed layer. The motion is there. Observations show that 
it dominates Eulerian measurements in the upper thermocline. The 
importance of this motion to mixed layer evolution needs to be 
established. 

C. SCOPE OF THIS STUDY 

The scope of this initial study shall be confined to a limited 
number of wave frequencies. There is insufficient time to concentrate 
on all the types of internal vertical motion, and all of the possible 
couplings and interactions in this thesis. The purpose of this initial 
research was to investigate the effect of dominant tidal-period internal 



15 



wave vertical motion on oceanic mixed layer dynamics, and to discern 
the effect upon the transient thermocline as observed during the mixed 
layer experiment. The effects of the vertical motion are explored 
using a one-dimensional bulk model of the mixed layer originally pro- 
posed by Garwood (1976,1977) and expanded here to include a prescribed, 
vertical velocity. 

The first hypothesis to be examined is can alteration of the 
existing mixed layer model be reasonably done? Can a relatively 
simple mathematical transformation be applied to include a prescribed 
vertical motion without altering the basic physics of the entrainment 
process? The study will show that this in fact can be done, as revealed 
by a careful examination of the interactions at each step of the alter- 
ation of the original model. 

A second hypothesis to be examined is that vertical motion produced 
by long period internal gravity waves for which the Rossby number is 
small, through interaction with cyclical boundary conditions (wind 
stress and heat flux) can produce two-dimensional (x-z) structures in 
the temperature field of a simple infinite ocean having horizontally 
homogeneous surface boundary conditions. An objective will be to 
demonstrate the frequency dependence of the internal vertical motion 
upon the phases of the cyclical heat fields, and the wind stress. It 
will be shown that there is a non-linear interaction between super- 
imposed diurnal and semi-diurnal internal tidal waves. Finally, this 
study will show the relative strengths of the resultant two-dimensional 
structures, and relate these to their importance in forming an improved 
thermodynamic ocean model. 
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II. MIXED LAYER MODEL AND VERTICAL MOTION THEORY 



A. MIXED LAYER MODEL 

To study the effects of internal vertical motion on mixed layer 
evolution, the Garwood (1976,1977) mixed layer model was used in a one- 
dimensional form. Garwood developed an ocean mixed layer model using 
the Navier-Stokes equation of motion with the geostrophic component 
eliminated, the continuity equation for incompressible water, the heat 
equation from the first law of thermodynamics, the conservation of salt 
equation, and a linearized equation of state. The entrainment hypothe- 
sis depends upon the relative distribution of turbulent energy between 
horizontal and vertical components and is offered as a plausible 
mechanism for governing both entrainment and layer retreat. 

Garwood suggested that planetary rotation influences the dissipa- 
tion of turbulence for deeper mixed layers and enables a cyclical steady 
state to exist on an annual basis. Furthermore, the rate of entrainment 
for the stable regime cannot be a simple linear extrapolation of the 
unstable situations. Unlike the atmospheric boundary case, most of the 
solar radiation does not penetrate the layer. Therefore, downward 
turbulent heat flux in the oceanic boundary layer is as important as 
the upward flux during the course of both diurnal and annual cycles. 

The non-linearity of the interface entrainment tendency parameter, 
which is greatest for stable surface boundary conditions, results in 
a modulation of the long term trend of the mixed layer depth by the 
diurnal component of surface heat flux. In this model, buoyant pro- 
duction is somewhat more efficient than shear production as a source 
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of energy for vertical mixing because of its unique effect on the 
vertical component of the turbulent velocity. 

"he governing heat equation for this model assumes two forms: 
one for the mixed layer, the other for below the mixed layer. For 
the mixed layer, the equation has the form 

d<T> = (w’T 1 ) - (w'T'L + qq /°Qdz , , 

~TT ^ p c p h '-h (1) 

where < > denotes the mixed layer average and Q is the net heat flux, 

a product of downward solar radiation, upward turbulent heat flux to 

the surface, and surface back radiation. Below the mixed layer the 

heat equation becomes 

dT = - d(wT) + ctqQ + K d 2 T 
d\ dz p C p dz 2 (2) 

where K is the eddy viscosity, and -3(wT)/iz is the term added for the 
vertical motion. 

The buoyancy flux at the entrainment zone is modelled according to 

asUTV = — <w’w‘>' ,2 <E> 

n h 

which is derived from the local turbulent kinetic energy (TKE) budget 
where 

z = -h (3a) 

In the entrainment zone, the buoyant damping attributed to (w T T r )_^ is 
assumed to be balanced by the convergence of TKE flux. This convergence 
is modeled as a function of the TKE kinetic energy components and the 
vertical distance over which this energy must be transported. <w 1 w ‘> 
and <E> are the bulk values of the vertical component and the total TKE, 
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respectively. The prognostic equations for these mixed-layer turbu- 
lence variables are derived using the bulk second-order closure methods 
of Garwood (1977), and are 

d (h<E> ) = - agh(w T T r L k + aghf(wT r )_ h + (VtY 1 

2 ' 2R* L 2 9J 

- (<E > ,/2 + fh) <E> ( 4 ) 



d (h<w'w'> ) = agh[(w'T') h + (w'T') ft ] + 

at 2 2 



where 



(<E> - 3<w r w I >)<E > l/2 - (<E> l/ 2 +fh)<E> (5) 

3 



u* = ( t //> 0 ) 1/2 

with r being the magnitude of the surface wind stress and p is the 

r Q 

air density, assumed to be a constant, and 

R* = aghAT 

(Au 2 + Av 2 ) 



is the bulk Richerdson number. The change in the total TKE, < r > , is a 
function of wind stress shear production, entrainment shear production, 
bouyant production/damping, and dissipation. The change in the vertical 
component, <w ‘ w *> , is a function of bouyant production/damping, pressure 
redistribution, and dissipation. The quantities AT, Au, and Av are the 
jump values between the mixed layer (assumed to be vertically uniform), 
and the underlying stable waters. These bulk TKE equations are solved 
algebraically assuming quasi-steady state for the TKE budgets 
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and 



_d_(h<E>) 

at 



0 



a (h<w'w'>) = 0 

at 



With this assumption the computed values of <w T w 1 > and <E> are in 
phase with the surface boundary conditions— the wind stress, ftu*, and 
the buoyancy flux, ag(w'T’ )g. 

Once the entrainment bouyancy flux is known from (3), the downward 
fluxes of momentum associated with entrainment at the base of the mixed 
layer, (w'u' )_ h and (w 1 v' are given by 

-(wV)_ h = w e Au = -(w'T'^Au ( 6 ) 

AT 

= » e Av “ -(w'T'LhJr (7) 



Here the entrainment velocity is 



AT 



and the entrainment velocity is also defined as 



— dh + w K 

e 17 " h 



(3) 

(9) 



if w a is greater than zero. If w g is not greater than zero, equations 
(10) and (11) are applied. 

1 . 3oundary Layer Deepening by Entrainment 

The entrainment heat flux (w'T 1 )_^ is determined algebraically 
from (3), and the steady-state forms of (4) and (5) in terms of given 
values of u*, (w'T 1 ) q , h and R*. This heat flux is then applied to 
the given temperature profile. The new h and new momentum and temperature 
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profiles are unique solutions provided that (i) the mixed layer is 
homogeneous, (ii) the value of h is just large enough to prevent 
unstable density profiles (dT/dz greater than zero), and (iii) heat 
is conserved. 

2. Boundary Layer Shallowing 

Shallowing occurs when there is adequate vertical turbulent 
energy to transport heat down to the base of the existing mixed layer. 
The depth at which the vertical heat flux (wT) vanishes establishes 
the new layer boundary, and will be the base of the new mixed layer. 

In this model, this occurs when (w'w 1 ) approaches zero. Then (3) is 
no longer applicable, and the steady state forms of (4) and (5) 
reduce to 

0 = mu? - (<E> ,/z +fh)<?> + qqh (wT\ 

* 2 ( 10 ) 

O = <E> VZ - (<E> ,/2 + fh)<E> + aqh (w'T T ) ft 

3 2 (11) 

Without the planetary rotation term, or if the downward surface buoyancy 
flux is dominant, h is proportional to the Obukhov length scale 

L =s u^ y /ag(w , T , ) 0 

If, on the other hand, (w'T' ) = 0, h is proportional to u*/f. In 

general, however, the depth of the shallowing mixed layer is a function 
of the two non-dimensional parameters, h/L and hf/u*. 

Heat and momentum are conserved in the profile specifications 
and, because the mechanical energy budget requires conservation of 
potential energy on the time scale of the turbulence, a numerical 
procedure is used that also preserves potential energy. This is 
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intended to ensure that the deepening rate is as correct as possible 
when the mixed layer deepens to its former position. 



B. VERTICAL MOTION THEORY 



By following a water particle as it is vertically advected by a 
prescribed velocity field, w(x,z,t), the distortion of any temperature 
profile, T(z), after a small, but finite length of time, (At), may be 
computed by non-linear theory. 

If we assume 

w = - Az co»(«t + kx) 

“F (12) 

is prescribed, then 

dz = - Az cos (cot + kx) 

dt 0 (13) 



This may be reformed and integrated 



= f -A cos(wt + kx) dt 

•*11 •'*,0 



where a water particle starting at time t-j at position z 1 is advected 
to position z 2 at time t 2 . Or following any constant isotherm 
initially at position Z] , to Z2 at t 2 gives 



z 2 = z, 



exp {_A_[ 3in(o/t 2 + kx) - sin(o>t, + kx)] } 
oj D 



(14) 



If t = t 2 - t ] « w* 1 , then 



z 2 = z, 



e xp { AAt[cos(o>(t| +J.) + kx)]J 



(15) 



or 



*2 S 




(16) 
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which is the general form approximated by the progression of equations 

in this study of internal vertical motion. Internal vertical motion 

is applied to the mixed layer model through equation (10). dh/dt, 

due to the internal vertical advection, w , , is derived from equation 

-h 

( 16 ). 

For future study, the salinity budget can be handled in a like 
manner. 
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III. EXPERIMENTS 



A. PRODUCTS 

The products used to identify features in the model runs and used 
for comparison are: 

1. Table (Sample Figure la) showing the hourly change in the mixed 
layer primary variables: 

H is the mixed layer depth in meters. 

T is the mixed layer temperature in degrees centigrade. 

S is the mixed layer salinity in parts per thousand. 

TAU is the surface wind stress in dynes per square centimeter. 

QNET is the net heat flux out of the mixed layer in watts 
per square meter. 

WDD is the value of the vertical motion parameter affecting 
the model in centimeters per second. 

DH is the change in mixed layer depth due to the vertical 
motion for that time period. 

2. Contours of the temperature variation, T(z,t), (Sample Figure 1) 

between z * 0 and z « -60 meters over the initial 24 hours of the model 

run. Varying resolution is used to better resolve the evaluation of 

T(z,t). Temperature contours are in degrees centigrade. 

3. Contours of temperature variation similar to 2, in a quasi two- 
dimensional mode. This is accomplished by running the model 24 times, 
each time passing the prescribed vertical motion through the 24 hour 
cycle of the model at a different, progressive phase. The 24 runs are 
then sorted by hour giving 24 data sets of 24 independent points with 

a spatial separation of kx from hour one through hour twenty-four. 
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DAY MR a 

1 L 6.3 
12 7.3 

1 3 8.0 

14 a. s 

15 9.9 

16 9.2 

17 9. S 

1 8 9.3 

1 3 10.1 

1 10 10.4 

1 11 10.7 

1 12 11.1 
1 13 11.4 

1 14 11.8 

1 IS 12.1 
1 16 12.4 

l 17 12'. 7 

1 18 13.0 

1 19 13.3 

1 20 13.6 

1 21 13.9 

l 22 14.1 

1 23 14.4 

1 24 14.7 

2 1 14.9 

2 2 IS. I 

2 3 IS. 2 

2 4 IS. 3 

2 5 IS. 3 

2 6 IS. 4 

2 7 iS.4 

2 8 1S.S 

2 9 15.6 

2 10 IS. 7 

2 11 IS. 9 

2 12 16.1 
2 13 16.3 

2 14 16. 5 

2 IS 16.6 
2 16 16.3 

2 17 17.0 

2 18 17.2 

2 19 17.4 

2 20 17. S 

2 21 17.7 

2 22 17.9 

2 23 18.1 

2 2 4 18.2 



T 3 

14.9S6 34.000 

14.939 34.000 

14.93S 34.000 

14.939 34.000 

14.949 34.000 

14.958 34.000 

14.969 34.000 

14.976 34.000 

14.978 34.000 

14.974 34.000 

14.963 34.000 

14.943 34.000 

14.926 34.000 

14.909 34.000 

14.893 34.000 

14.877 34.000 

14.362 34.000 

14.847 34.000 

14.833 34.000 

14.819 34.000 

14.306 34.000 

14.793 34.000 

14.781 34.000 

14.767 34.000 

14.762 34.000 

14.762 34.000 

14.768 34.000 

14.777 34.000 

14.789 34.000 

14.901 34.000 

14.914 34.000 

14.823 34.000 

14.829 34.000 

14.830 34.000 

14.825 34.000 

14.814 34.000 

14.803 34.000 

14.793 34.000 

14.793 34.000 

14.773 34.000 

14.763 34.000 

14.753 34.000 

14.743 34.000 

14.734 34.000 

1 4 . 72 S 34.000 
14.716 34.000 

14.707 34.000 

14.697 34.000 



TAU £NET 

0. SO - 3. S 

O.SO -100.8 
O.SO -184.3 
O.SO -248.4 
0.50 -288.7 

O.SO -302. S 
O.SO -288.9 
O.SO -248.7 
0.50 -184.7 

O.SO -101.3 
O.SO - 4.1 

O.SO 100.2 
O.SO 100.8. 
O.SO 100.8 
O.SO 100.8 
O.SO 100.8 
O.SO 100.8 
O.SO 100.8 
O.SO 100.3 
O.SO 100.8 
0.50 100.8 

O.SO 100.8 

0.50 100.8 

O.SO 100.8 

O.SO - 3.S 

0.50 -100.3 

O.SO -184.3 
O.SO -248.4 
O.SO -288.7 
0.50 -302. S 

O.SO -288.9 
O.SO -248.7 
0.50 -184.7 

O.SO -101.3 
O.SO - 4.1 

O.SO 100.2 
O.SO 100.3 
0.50 100.8 

O.SO 100. a 
O.SO 100.8 
O.SO 100.8 
O.SO 100.8 
O.SO 100.8 
O.SO 100.8 
0.50 100.8 

O.SO 100.8 
O.SO 100.8 
0.50 100.8 



MOD DB 

0.0 0.0 

0.0 0.0 

0.0 0.0 

0.0 0.0 

0.0 0.0 

0.0 0.0 

0.0 0.0 

o.o o.o 

0.0 0.0 

0.0 0.0 

0.0 0.0 

0.0 0.0 

0.0 0.0 

0.0 0.0 

0.0 0.0 

0.0 0.0 

0.0 0.0 

0.0 0.0 

0.0 0.0 

0.0 0.0 

0.0 o.o 

0.0 0.0 

0.0 0.0 

0.0 0.0 

0.0 0.0 

0.0 0.0 

0.0 0.0 

0.0 0.0 

0.0 0.0 

0.0 0.0 

0.0 0.0 

0.0 0.0 

0.0 0.0 

0.0 0.0 

0.0 0.0 

0.0 0.0 

0.0 0.0 

0.0 0.0 

0.0 0.0 

0.0 0.0 

0.0 0.0 

o.o o.o 

0.0 0.0 

0.0 0.0 

0.0 0.0 

0.0 0.0 

0.0 0.0 

0.0 0.0 
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6.2613 
7.2840 
9.0111 
9.4966 
9.9230 
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9.S392 
9.8356 
10. 1260 
10.4177 
10.7468 
11. 1046 
11.4344 
11.7720 
12. 1007 
.2.4024 
2.7107 
3. 3222 
13.2981 
13.5796 
13.5668 
14. 1371 
14.3931 
14.6S36 




Figure 1. (a) Table of hourly values of 

problem parameters. H is the entrainment 
depth in m. T is mixed layer temperature 
in °C. S is mixed layer salinity in 0 / 00.7 
TAU is the surface wind stress in dy^es/cm". 

Q.NET is the net heat flux in watts/nM WDD 
is the vertical motion parameter. DH is the 
hourly change in mixed layer depth due to 
the vertical motion in cm. (b) The contours 
of T9z,t) as evolved over a 24 hour period to 
a depth of 60 meters. (c) Successive temperature 
profiles from hour one through hour twenty four. 
The first profile is true to the temperature 
scale on the graph. Succeeding profiles are 
offset. 
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It is assumed that there is a two-dimensional internal wave in x and 
z passing through this ocean system. Heat flux and wind stress are 
taken to be the same at all points, with no interaction between any 
of the points (advection is limited to the vertical). 

4. Progressive temperature profiles, T(z), to a depth of 60 meters 
(Sample Figure 1c). Only the first profile matches the horizontal 
temperature scale which is degrees centigrade. The remaining profiles 
have been consecutively offset. The accompanying table has the mixed 
layer temperatures in degrees centigrade, and the mixed layer depth 

in meters. 

5. Progressive temperature profiles as in 4, but in the two- 
dimensional mode of 3. Profiles are offset from the first, again, to 
provide an insight to the two-dimensional nature to the data. 

B. STANDARDS 

Prior to the examination of the addition of vertical motion to the 
mixed layer model (MLM), it is beneficial, here, to examine the results 
of running the mixed layer model in the two modes which will be used 
as the standards for comparison with the vertical motion modelling 
experiments. Both standard cases are runs of the mixed layer model 
from which the effects of evaporation and precipitation have been 
removed. This has been done to aid in the interpretation of later 
results by limiting the number of interaction variables to a manageable 
quantity. The first standard is a run of the MLM in the primarily wind- 
driven deepening mode. The second standard is for the MLM in a cycli- 
cal shallowing/deepening response to a strong downward heat flux. 
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1 . Standard 1 ($1 ) 



This run was initialized with a temperature profile having a 
surface and mixed layer temperature of 15.0°C. At 5.0 m the initial 
temperature decreases linearly with the form 

T(z) = T, - 5.0 X I0" 4 (Z- H) } 

Where Z is the absolute value of z in centimeters, and H is the mixed 
layer depth in centimeters. This profile extends to 200m. The result- 
ing profile is shown in Figure 2. 

The surface wind stress (TAU) is a constant 0.50 dynes per 
square centimeter for the entire run. Net heating reaches a maximum 
of 302.5 watts per square meter at hour six (local noon). Net heat 
loss after sunset is a constant 100.8 watts per square meter. 

With no vertical motion the layer deepens continuously over 
the two-day period of the model run (Figure 1), as expected for pri- 
marily wind-driven mixing. In examining the first twenty-four hour 
period it is seen that the mixed layer deepens from an initial depth 
of 4.5 m to a depth of 14.7 m. The layer cools from 15°C to 14.767°C 
over the period, with a short period of mixed layer temperature 
increase of six hours related to that time when the net heat flux 
is sufficient to overcome the cooling attributed to the entrainment 
of cooler water into the layer. 

The salinity in the layer remains constant because there is 
no precipitation or evaporation allowed which could change the surface 
salinity, and the below layer salinity is assumed to be a constant 
34.00 part per thousand. Therefore, there are no sources of water 
with a salinity other than 34.00 parts per thousand. 
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Figure 2. Initial temperature and salinity profiles for Case SI. 
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The downward solar radiation, Q , increases from sunrise (hour 1) 
to noon (hour 6), and then decreases to zero at sunset (hour 12.) A 
constant is assumed for the net upward turbulent heat flux to the atmos- 
phere, plus back radiation of 100.8 watts per square meter for the 
entire day. 

The one-day response to this cycle can be seen in the contours 
(Figure lb). The nixed layer deepens continuously over the first 24 
hours. 3elow the layer temperature, at any given depth not penetrated 
as the layer deepens, remains constant over time. In the layer the 
temperature is seen to cool slightly initially, then warm to a maximum 
value at about hour nine (1500 local time), and then cool again for 
the remainder of the period. 

In the T(z) graph this sequence is seen again (Figure 1c). 

Here it is seen that all of the slopes of the profiles below the layer 
are the same from hour one to twenty-four, indicating that with no 
vertical motion, the characteristics below the layer will not change 
with time. The included table shows the relationship of the temper- 
ature, in degrees centigrade, and the layer depth, in meters, and 
indicates the rate and magnitude of the deepening. 

The rate of deepening (Table 1) reflects the surface boundary 
conditions at any given time as well as layer depth and stratification 
below. As the run begins, the layer deepens quickly in response to 
the surface wind stress and the short distance over which it must act, 
as seen in equations (3) and (4), At hour one the layer deepens 
beginning at a shallow value where the wind flux has a high impact 
from equation (4), during a period where the heat flux into the layer 
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is a minimum. As the heat flux increases and the layer deepens the 
wind stress is relatively less effective in deepening, and the entrain- 
ment now effected by equation (5) slows until hour nine when the solar 
radiation goes to zero, and the deepening rate increases. This increase 
is relatively small because of the relatively large temperature jump 
(AT) in the entrainment zone by this time. After hour eleven (1800 
local time) the minor heat loss, small wind stress and increased depth 
involved all act to slowly decrease the deepening rate. 



Table 1. Deepening Rate in Meters Per Hour for Case SI 



Hour 


Deepening 

Rate 


Hour 


Deepening 

Rate 


i 


1.7613 


13 


0.3296 


2 


1.0227 


14 


0.3376 


3 


0.7271 


15 


0.3287 


4 


0.4877 


16 


0.3017 


5 


0.4242 


17 


0.3083 


6 


0.3221 


18 


0.3115 


7 ' 


0.2941 


19 


0.2759 


8 


0.2964 


20 


0.2815 


9 


0.2904 


21 


0.2872 


10 


0.2917 


22 


0.2703 


11 


0.3291 


23 


0.2560 


12 


0.3580 


24 


0.2607 



2. Standard 2 (S2 ) 

This run is begun with the same initial layer temperature as 
for standard 1. Here, however, the initial mixed layer depth is taken 
to be 49,5 meters. At 50.0 meters, the profile again assumes the 
linear variability of equation (18), extending to 200.0 meters 
(Figure 3). 

All other parameters are the same as those of SI except for 
the prescribed surface heat flux. The net heat flux has increased to 
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Figure 3. Initial temperature and salinity profiles for Case S2. 



a maximum of 504.2 watts per square meter by assuming an extremely 
strong solar heating maximum (Q^) of 600 watts per square meter at 
local noon. 

In this case the layer depth remains at 49.5 meters for the 
first hour (Figure 21a) because the wind stress is not strong enough, 
coupled with the available heat energy, to cause further entrainment 
at this depth, and the solar heating is not yet strong enough to cause 
(w'w* ) to go to zero so that equations (10) and (11) are applied. By 
the second hour the net downward surface heat flux and the surface 
wind stress are sufficient to cause (w'w* ) to go to zero, and a new h 
is formulated through the solution of equations (10) and (11), forming 
a new surface mixed layer with a depth of about 18 m. The layer then 
continues to shallow until hour 6 (local noon) in response to the 
increased downward heat flux. Stirring by the wind stress is relatively 
small as its effect is overcome by the buoyancy flux associated with 
net surface heat flux. The layer temperature increases as long as 
there is heating sufficient to overcome the minor cooling due to the 
entrainment of cooler below layer water which resumes at local noon, 
but is not strong until hour 10 (1600 local time). At that time the 
layer cools due to the combined effects of the surface heat loss, and 
the entrainment of the cooler below layer waters. 

The salinity, again, remains constant because there is no 
source of water of a salinity other than 34.00 parts per thousand. 

As shown by the contours (Figure lib) the original layer of 50 meters 
shallows to hour 6, not forming a distinct layer, since equations 
(10) and (11) are applied at each hourly step, until hour 7. The 
transient thermocline prior to hour 7 is weak. This is seen more 
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clearly, perhaps, on the T(z) profiles (Figure 11c). The originally 
sharp interface at 50 meters becomes diffuse with time due to small 
(K-j. = 0.2 cm^/s) below-layer diffusion. This is clearly seen in the 
T(z) profiles, and it is the reason for the non-horizontal contours 
near the 50 m depth. The variability of the temperature structure 
below the newly formed surface mixed layer is a function of the shal- 
lowing adjustment by equations (10) and (11) as applied from hour two 
through six. 

Again it is evident in the T(z) profiles that the profiles 
keep the same slope below the layer except for the originally sharp 
interface at 50 meters where diffusion of the layer occurs with time. 
This diffusion area is also seen to hold at a relatively constant 
depth with time. 

The T(z,t) evolution reflects the turbulent kinetic energy 
budget for this problem (Table 2). The initial decrease in h from a 
deep layer to a shallow layer is clearly the result of the surface 
heating and the surface wind stress acting to form a new, warmer, 
shallower surface layer. Energy is increased in the surface mixed 
layer due to excessive solar insolation during the day. The Obukhov 
length is reduced from hour two through six, and the newly forming 
mixed layer adjusts accordingly. By hour seven the solar insolation 
begins to reduce. As this happens, the Obukhov length increases, 
allowing the newly formed layer to deepen. Deepening does not keep 
up with the Obukhov length because of the damping due to the energy 
transport which must take place during the mixing process. By hour 
11 the net heat flux is negligible decreasing to a constant surface 
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loss after sunset. The relative magnitude of the depth increase 
(dh/dt) with time when compared with that of SI is due to the strength 
of the temperature jump (AT) at the layer, and the available energy 
in the layer. In SI a AT of the order of ,4°C was encountered, 
while in S2, due to the slow formation of the new surface layer, a 
AT of the order ,08°C is the maximum encountered. This allows for 
more rapid deepening. The rate of deepening increases after hour 16 
(2200 local time) because, as seen in the T(z) profiles (Figure lie) 
the AT at the entrainment zone is decreasing with time due to the net 
upward surface heat flux. 

Table 2. Deepening Rate in Meters Per Hour for Case S2. 



Hour 


Deepening 

Rate 


Hour 


Deepening 

Rate 


1 


0.0152 


13 


0.8464 


2 


-31.4079 


14 


0.8317 


3 


-1.9247 


15 


0.8252 


4 


-3.1646 


16 


0.8235 


5 


-0.9575 


17 


0.8278 


6 


-0.2935 


18 


0.8404 


7 


0.3766 


19 


0.8596 


3 


0.3445 


20 


0.8879 


9 


0.4944 


21 


0.9186 


10 


0.5383 


22 


0.9516 


11 


0.6793 


23 


0.9370 


12 


0.8371 


24 


1.0256 



C. EXPERIMENTS 

The addition here to the existing theory of Garwood (1976-77) is 
approached in a step-by-step procedure, gradually increasing the com' 
plexity of the problem. A quick reference for this progression can 
be seen in Table 3. A basic assumption that was made in all of the 
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able 3. Ini tial conditions, boundary conditions, and progression 
of vertical motion inputs for all experiments. 



For all Cases: 

1. Net surface heat loss is constant at 100.8 watts/m 2 . 

2. Surface wind stress is constant at 0.S dynes/cm 2 . 

3. The initial layer temperature is 15.0*C. 

4. The salinity is constant to 200 m at 34.00°/oo. 

5. No precipitation or evaporation is allowed during the run. 

6. The vertical motion is exponentially decreasing to the surface. Maximum motion 
is held at 36.0 cm/hr at 10 m, except for Case IVb, where it is 18.0 cm/hr at 
10 a. 



Case 


B. C. 

Qs aax. 
(watts/m*) 


I. 

Temp. 

Profile 


C. 

MLD 

Cm) 


Vertical Motion 


SI 


400.00 


Fig. 2 


4. S 


None 


S2 


600.00 


Fig. 3 


49.5 


’ None 


la 


SI 


SI 


SI 


Constant downwelling during model run. 


lb 


S2 


S2 


S2 




Ic 


SI 


SI 


SI 


Constant upwelling during model run. 


Id 


S2 


S2 


S2 




Ila 


SI 


SI 


SI 


Pulsed downwelling, 12 hour cycle. Phase 
shifted to peak at hours 5, 18, and 24. 


lib 


S2 


S2 


S2 


lie 


SI 


SI 


SI 


Pulsed upwelling, 12 hour cycle. Phase 
shifted to peak at hours 5, 18, and 24. 


lid 


S2 


S2 


S2 


Ilia 


SI 


SI 


SI 


Single diurnal wave passage. Phase shifted 
to give maximum downward motion at all hours. 


mb 


S2 


S2 


S2 


Hours 1, 7, 13, and 19 discussed. 


IVa 


S2 


S2 


S2 


Continuous diurnal wave system passage. Same 
phase shifts as Case III. Hours 6, 12, 18, 
and 24 discussed. 


IVb 


S2 


S2 


S2 


Continuous semi-diurnal wave system passage. 
Phase shifted as in Case III. Runs with peak 
motions at hours 6 and 18, and 12 and 24 are 
discussed. 


Va 


S2 


S2 


S2 


Superimposed diurnal and semi-diurnal continuous 
wave system passage. Both waves in phase. 

Input phase shifted as in Case III. 


Vb 


S2 


S2 


S2 


As in Case Va, but phases of diurnal and semi- 



diurnal waves varied relative to one another. 
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work was that the vertical motion moved the water up and down, but did 
not directly induce turbulence (i.e., mixing). This is a valid 
assumption if no shear is associated with this motion. 

The cases of progressively complicated vertical motion examined 
were: 

1. Linearly (in z) varying vertical motion, constant over the 
entire period. 

a. Upwelling 

b. Downwelling 

2. Linearly varying vertical motion pulsed in time, 
a. Short daily upwel 1 ing/downwelling events. 

3. Linear internal wave motion. 

a. Passage of a single wave pulse through the experimental 
point over the 24 hour daily cycle. 

4. Non-linear wave variations coupled with diurnal and semi- 
diurnal boundary conditions (heating, wind). 

5. Non-linear response to superimposed waves of different periods. 
Diurnal and semi-diurnal tidal waves are combined and phase- 
shifted to determine non-linear interactions. 

1 . Case I 

Both upwelling and downwelling are realistically induced with 
the realization that in a mixed homogeneous layer the vertical dis- 
placement due to a linear w(z) can be shown to be exponential, going 
to zero at the surface. Mathematically we know that 



= w(z) 
dt 



w(z) = Az 
D 

where D is a scaling depth, it follows that 



( 13 ) 
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dz 

dt 



Az 

0 



or 




Integration from t^ to t 2 yields 



z, e 




(19) 



which is the equation upon which the following work is based, and which 
agrees with equation (16) derived earlier. 

The temperature and salinity profiles must reflect "correct" 
variation if equations (16) and (19) are correct. Upwelling and/or 
downwelling will not cause isolines to break the surface or each other. 
Mixing must be done to induce this. Therefore, above the mixed layer 
depth T new and $ new must equal T^ and when there is vertical 
motion without mixing or surface exchange. Below the entrainment depth 
the profiles must be altered to reflect the type of vertical motion 
(upward or downward) in accordance with equation (16). Figure 4 shows 
the expected result from upwelling on a temperature profile. The mixed 
layer depth has decreased from time one to time two (after upwelling), 
and the slope of the profile below the mixed layer has been decreased. 
In all the cases the AT at the entrainment zone has been preserved. 

a. Constant Downwelling (Case la) (Figure 12) 



(See Table 3). A constant downward vertical motion is introduced giving 
the vertical motion parameter, WDD, (WDD = A t/D) a value of .03600, 
which equates to a value of vertical motion of .01 cm/s or 36 cm/hr at 



The initial conditions used are the same as those for SI 
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Figure 4. Temperature profile variation from tj to t 2 under upwelling conditions. 



ten meters. This is a rather large value, but it demonstrates the inter- 
actions well. The vertical motion is applied beginning the first hour, 
and is not applied as the initial condition on the temperature and 
salinity profiles. 

The result of this downward motion is, as expected, a more 
rapid deepening of the surface mixed layer (Figure 12a) when compared 
to SI. Comparing the depth change due to the mixing in this case (Table 
4) to the SI mixing, the same sequence of events as described for SI is 
seen, and the values are comparable, indicating that the physics of the 
model has not been altered by the addition of the vertical motion para- 
meter in the downwelling mode. It is noted that the effect of the 
vertical motion quickly becomes greater than the effect of mixing due 
to the large value of vertical mixing assigned in this run of the model. 



Table 4. Deepening Rate in Meters Per Hour for the 
Mixing aad Vertical Motion of Case la. 



Hour 


Mixing 

Deepening 

Rate 


Vertical Motion 
Deepening 
Rate 


Hour 


Mixing 

Deepening 

Rate 


Verticle Motion 
Deepening 
Rate 


1 


1.7613 


0.2295 


13 


0.3201 


0.5729 


2 


0.9705 


0.2735 


14 


0.3169 


0.6055 


3 


0.6606 


0.3077 


15 


0.3134 


0.6392 


4 


0.4759 


0.3365 


16 


0.3102 


0.6740 


5 


0.3535 


0.3613 


17 


0.3068 


0.7099 


6 


0.2945 


0.3858 


18 


0.3034 


0.7471 


7 


0.2444 


0.4089 


19 


0.3017 


0.7855 


3 


0.2282 


0.4323 


20 


0.2983 


0.8252 


9 


0.2394 


0.4569 


21 


0.2949 


0.8663 


10 


0.2513 


0.4828 


22 


0.2914 


0.9087 


11 


0.2793 


0.5108 


23 


0.2881 


0.9526 


12 


0.3219 


0.5413 


24 


0.2094 


0.9952 



The sequence is vividly displayed in the contours (Figure 12b) 



where it is seen that the temperature contours deepen as time progresses 
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as one would expect in a downwelling case. The surface heating due 
to the incoming solar radiation is noted again, followed by cooling 
as the sun sets and heat is then lost at the surface and colder water 
is entrained at the bottom of the mixed layer. It is also noted that 
in only twenty-four hours, the vertical motion added here has deepened 
the layer to almost twice the value of SI. 

The T(z) plot (Figure 12c) demonstrates the deepening 
involved and the effects on the temperature profile below the mixed 
layer. Note the increased slope (dT/dz) with time, as warmer near- 
surface temperatures (water) are carried downward by the downwelling. 
Also, one can see the effect of the preservation of the AT profile at 
the entrainment zone from the mixing step through the vertical motion 
progression in the model where at times there is an apparent discon- 
tinuity in the profile. This does not appear to cause a problem in 
the execution of the model, and will be seen in many of the following 
T(z) plots. One can further see that the vertical motion addition has 
caused a slight decrease in the jump strength of AT at the layer over 
that in SI. This can be attributed to the reduced mixing for the 
larger mixed layer depths. Where, with no downwelling a AT of 1.0°C 
may have occurred in 18 meters, it now occurs at a greater depth. 

This effect is transmitted into the entrainment zone as this below 
layer water is mixed in, causing the AT to be reduced. 

d. Constant Downwelling (Case lb) (Figure 13) 

The initial conditions are those of S2 (See Table 3). A 
constant downward motion is introduced equal to that used in the 
previous case. 
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The result of the downward motion is again, as expected, 
a more rapid deepening of the layer after the new surface mixed layer 
is formed through the surface wind stress and heat flux as applied by 
equations (10) and (11) (Figure 13a). Comparing the depth change due 
to the mixing here (Table 5) to the S2 mixing shows that the values 
are of the same order, which shows that the energy relations of the 
original model are not impaired by the addition of vertical motion 
in this mode of the model either. Here it is noted that mixing drives 
the production of a new surface mixed layer beginning at hour 2 (0800 
local) (Figure 13b). The shallowing due to the surface heating is 
dominant until hour seven (1300 local) when the downwelling becomes 
the dominant factor until hour nine. From hour nine to hour fifteen 
the effects of mixing and downwelling are of the same order. After hour 



Table 5. Deepening rate in meters per hour for the 
Mixing and Vertical Motion in Case lb. 



Hour 


Mixing 

Deepening 

Rate 


Vertical Motion 
Deepening 
Rate 


Hour 


Mixing 

Deepening 

Rate 


Vertical Motion 
Deepening 
Rate 


1 


0.0152 


1.8150 


13 


0.8093 


0.6986 


2 


-33.4052 


0.6571 


14 


0.7343 


0.7511 


3 


-2.5282 


0.5885 


15 


0.7445 


0.8059 


4 


-3.7245 


0.4735 


16 


0.7131 


0.8616 


5 


-1.4432 


0.4392 


17 


0.7029 


0.9190 


6 


- .7286 


0.4286 


18 


0.6768 


0.9775 


7 


0.1233 


0.4489 


19 


0.6716 


1.0379 


8 


0.3275 


0.4773 


20 


0.6628 


1.1002 


9 


0.4266 


0.5104 


21 


0.6420 


1.1641 


10 


0.5787 


0.5504 


22 


0.6316. 


1.2299 


11 


0.6446 


0.5942 


23 


0.6224 


1.2978 


12 


0.8003 


0.6986 


24 


0.6120 


1.3678 



fifteen the downwelling again, slowly, becomes the dominant factor for 
layer deepening. 



It is seen that the layer temperature, again, immediately 
begins to warm (Figure 13c), and in fact, while the dominant factor is 
the large amount of surface heating, the layer warms at exactly the 
same rate as in S2, After hour ten the layer again begins to cool, but 
because of the deepening of the mixed layer due to downwelling, the 
mixing of the surface wind stress is reduced. This, along with the 
fact that a deeper layer has available more mass of water to dilute 
the effect of the entrained water, results in less cooling over the 
remaining twelve hours. 

Because of the downwelling effect on diffusing the temper- 
ature gradient it is also seen (Figure 13b) that the entrainment zone 
(layer boundary) remains rather diffuse until hour eight, while in S2 
it was sharpened by hour six. This helps explain the magnitude of the 
deepening effect in the period from hour nine to fifteen, though the 
depth of the entrainment zone is increasing faster than that in S2 and 
reaches a difference of five meters from that in S2 by hour fifteen. 

The AT across the bottom of the layer for the downwelling case is 
slightly less than that in S2, allowing the mixing to be a more impor- 
tant factor for deepening in this case. The following decrease in the 
mixing is then due to the ever increasing depth through which the mix- 
ing must take place. The wind stress must act over distances that are 
twice as large as those in S2, and although the AT at the layer depth 
is not as great, the effect of the depth is quite noticeable. 

The original deep layer is seen to begin to diffuse, but 
is quickly carried below the 60 meter level of the graphs. As stated 
above, the new surface mixed layer begins to show a significant AT 



42 



at hour eight. Extremely close examination of the T(z) profiles for S2 
and the profiles for this run does show a smaller AT at the entrainment 
zone for this case during the period from hour seven to hour sixteen. 

After hour 16 there is no discernable difference in the AT in the profiles, 
c. Constant Upwelling (Case Ic) (Figure 14) 

The initial conditions are those of SI . A constant upward 
vertical motion is introduced beginning at hour one giving the vertical 
motion parameter, WDD, a value of .03600, or about 36 cm/hr at ten meters. 
Again, this is high, but it is done with the realization that later work 
will be with wave motion where this will not be an unrealistic value, 
and using it here provides a good benchmark for the effects produced by 
this motion. 



Table 5. Deepening Rate in Meters Per Hour for the 
Mixing and Vertical motion of Case Ic. 



Hour 


Mixing 

Deepening 

Rate 


vertical Motion 
Deepening 
Rate 


Hour 


Mixing 

Deepening 

Rate 


Vertical Motion 
Deepening 
Rate 


1 


1.7613 


-.2214 


13 


0.4239 


-.3026 


2 


1.1031 


-.2526 


14 


0.4528 


-.3079 


3 


0.6746 


-.2675 


15 


0.4845 


-.3142 


4 


0.6103 


-.2796 


16 


0.4774 


-.3199 


5 


0.4954 


-.2873 


17 


0.2909 


-.3189 


6 


0.3417 


-.2892 


18 


0.2871 


-.3178 


7 


0.3184 


-.2902 


19 


0.3046 


-.3173 


8 


0.3097 


-.2909 


20 


0.3223 


-.3175 


9 


0.3133 


-.2917 


21 


0.3407 


-.3183 


10 


0.3254 


-.2929 


22 


0.2596 


-.3163 


11 


0.3435 


-.2947 


23 


0.2739 


-.3148 


12 


0.3938 


-.2982 


24 


0.2891 


-.3138 



The result of the addition of this upward motion is, as 
anticipated, a slowing of the deepening rate of the mixed layer (Figure 
14a) when compared to SI. In fact, due to the size of the vertical 
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motion parameter, an effect is produced which almost holds the mixed 
layer at a constant depth. This is displayed more dramatically after 
two model days have been executed, and can be seen to reach a near 
equilibrium condition after four to five days. Comparing the deepening 
rates due to the effect of the mixing in the model and the vertical 
motion in this case (Table 5) to the 51 mixing it can be seen quite 
clearly that the effects described for the mixed layer model standard, 
SI, are still present. From equations (1) and (3) the slightly 
increased magnitude of the cooling of the mixed layer by entrainment 
and the subsequent heating are attributed to the reduced depth over 
which the wind stress and heat flux must be distributed. Again, it 
can be seen, that the model TKE balance, equations (3-5), is not 
violated by the addition of this upwelling motion. 

The layer temperature initially falls (Figure 14a), then 
rises until 1400 local time. It then falls again in response to the 
mixing and the surface heating cycle, as in SI. The heating of the 
layer is not as great here as in the downwelling case (Figure 13), 
because the entrainment is sufficiently enlarged due to the smaller h, 
and because cooler water becomes available for entrainment more quickly 
due to the upwelling. This is more clearly seen in comparing the 
contoured fields, T(z,t), (Figure 14b) of this case with those of the 
downwelling case (Figure 13b). By the end of the twenty-four hour 
period this has resulted in a very strong temperature jump, AT, at 
the entrainment zone, as is seen best in the T(z) profiles (Figure 14c). 
It is seen that AT grows with time. 

The strengthening of the jump, and the lifting of the 
cooler water is demonstrated clearly by the contours, and it is 
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clearly evident the contours not only lift, but also merge in time, 
indicating the AT grows with time. This effect is seen in the T(z) 
profiles by the increased gradient of the profiles with time. 

Equilibrium is reached in this case as the jump gets too 
large for the mixing to have an effect greater than the effect of 
the vertical motion. Equation (8) shows that dh/dt vanishes as w e 
approaches w.^. 

d. Constant Upwelling (Case Id) (Figure 15) 

The initial conditions applied are identical to those of 
S2. A constant upward vertical motion is introduced at hour one of 
about 36 cm/hr at ten meters, as before. 

The net result of the vertical motion after 24 hours is a 
reduced layer depth in comparison to the case having no vertical motion 
(S2). The total reduction in depth of six meters, however, is not as 
great as the previous case (Ic). After two days it is seen that the 
mixed layer is approaching an equilibrium depth of near eleven meters 
(Figure 15a). 

The deepening rates (Table 6), when compared to those of 
S2, demonstrate that the basic model is still preserved in this run. 
From hours two through five the model is in the shallowing mode, and 
the new layer depth is established by equations (10) and (11), while 
there is a massive influx of solar radiation to the sea surface. The 
increased deepening rate due to the mixing seen from hour twelve 
through twenty-three is attributed to the effect of equations (1), 

(4) and (5). The mixed layer is maintained at a shallow level by the 
upward vertical motion, therefore, the mixing takes place over a much 
smaller distance than in S2. 
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Table 6. Deepening Rate in Meters Per Hour for the 
Mixing and Vertical Motion of Case Id. 



Hour 


Mixing 

Deepening 

Rate 


Vertical Motion 
Deepening 
Rate 


Hour 


Mixing 

Deepening 

Rate 


Vertical Motion 
Deepening 
Rate 


1 


0.0152 


-1,7509 


13 


0.9465 


-.4751 


2 


-29.4675 


-0.6470 


14 


0.9470 


-.4917 


3 


-1.3404 


-0.5767 


15 


1.0159 


-.5103 


4 


-2.6147 


-0.4639 


16 


1.0611 


-.5298 


5 


-0.5166 


-0.4292 


17 


1.1449 


-.5515 


6 


0.1312 


-0.4186 


18 


1.2301 


-.5755 


7 


0.3085 


-0.4148 


19 


1.3043 


-.6013 


8 


0.4578 


-0.4163 


20 


1.4087 


-.6298 


9 


0.5885 


-0.4224 


21 


1.5577 


-.6626 


10 


0.6775 


-0.4314 


22 


1.6638 


-.6980 


11 


0.7113 


-0.4413 


23 


1.3743 


-.7219 


12 


0.9078 


-0.4578 


24 


0.6909 


-.7208 



Further, there is a concentration of energy in this shallower 
layer. The great variance in the last hour is a direct result of the 
magnitude of the upwelling, which, by hour four, has set up a large 
gradient in the temperature (Figure 15b). By hour twenty-four a gradient 
has been set up with such a large slope that mixing is becoming very 
difficult. Continued upward motion of this magnitude, without shallow- 
ing as is seen on the second day, would quickly make it extremely 
difficult to entrain much water into the mixed layer. 

The progression in the change of the mixed layer temperature 
(Figure 15a) shows that as the model is shallowing, governed by equation 
(10) and (11), the same adjustment takes place as is seen in S2. The 
temperatures produced are exactly those seen in S2. This is another 
indication that the model is performing as before. Since w e is zero, 
the upwelling has no effect on the setup of the new layer. If the H 
values had been listed immediately after the mixing phase of the compu- 
tation, this could be seen more clearly. After the layer again begins 
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to mix downward, however, the upwelling begins to have some effect on 
T and h due to the upward transport of the entrainment zone. Though 
the layer is one to two meters shallower than in S2, the temperature 
of the layer is the same as that in the deeper standard to hour ten. 

At hour ten the heat flux of equation (1) has been reduced sufficiently 
so that wind mixing again becomes the dominant deepening force. At 
this point the upwelling has moved water cooler than that in S2 at the 
same time to the entrainment zone. The change in entrainment rate due 
to the upwelling, the entrainment of cooler waters, and the net surface 
heat loss after hour eleven combine to produce a cooling of ,147°C from 
the maximum warm temperature, as compared to ,109°C for S2. For the 
mass of water involved, and the relatively short time involved, this 
is a significant difference. 

The effect of the upward vertical motion is perhaps best 
shown in the T^z) diagram (Figure 15c), where it is seen that the model 
initially adjusts, through equations (10) and (11), to a new surface 
layer which slowly strengthens, in terms of the jump, AT, as it en- 
trains. As the process continues the profiles below the layer rise 
by upwelling and the initial deep interface is smoothed by diffusion 
until it disappears. By hour twenty- four a sufficiently stable pro- 
file has been set up to make deepening of the layer more difficult. 

This indeed is the case in the second day of the run (Figure 15a) as 
the layer reaches a temporary equilibrium. It is temporary, in that 
once the profile has been reached which makes deepening sufficiently 
difficult through equations (4) and (5), the dominant force becomes 
the upwelling, w_^ of equation (9). 
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The contours (Figure 15b) further demonstrate the effect 
of the upward lifting of the cooler water, as the contours tend to 
converge with time. The initial interface is quickly lost, as the 
model smooths this into its uplifted profile. 

2. Case II 

The next step of this progressive series of tests is to explore 
the effects of a simple pulsed motion varying in time. This is accom- 
plished by varying the input motion in a sinusoidal pattern. By using 
the equation 

w ( t ) = AA sin ( 2irt + kx) 

° . ■ (20) 

as the particular form of equation (12) applied for this case, where 
AA is the velocity scale factor such that AA = Az/D is scaled for an 
appropriate motion at 100 meters, and limiting the value of w Q (t) to 
only values in one direction (+ or -), and zero, a pattern of motion 
can be produced such as that seen in Figure 5. A pulse of vertical 
motion of twelve hour duration is produced, followed by a twelve hour 
period of no motion. To determine the relationship of this motion 
with the cyclical boundary conditions one need only shift the phase 
of the pulse relation. This will cause the motion to peak at varying 
selected times during any twenty-four hour period run. 

a. Pulsed Downwelling (Case 1 1 a) (Figure 16, 17 and 13) 

The initial conditions are the same as those of SI. Three 
different phase relationships of the downwelling to the cyclical 
boundary conditions are presented. In the first (Figure 16) the 
maximum of the downwelling occurs 180° out of phase with the heating 
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Figure 5. Pulsed upwelling form (in relation to model time/hrs at point x. Amplitude is variable 



cycle (that period of the model defined as night, or the cooling cycle). 
In the second (Figure 17), the downwelling occurs during the last seven 
hours of night, and during the first five hours of daylight. The max- 
imum downwelling occurs six hours before the heating maximum. Finally, 
in the third case (Figure 18), the downwelling occurs in phase with 
the maximum heating cycle. In all three cases, the downwelling event 
is a twelve hour pulse with a maximum value of 36 cm/hr at ten meters. 

(1) 180° Out of Phase (Figure 16) . For this case it is 
seen that the model evolves exactly as does SI until hour thirteen 
when downwelling begins (Figure 16a). All the fields produced are 
exactly the same as in SI until hour thirteen for which the evolution 
has been discussed. 

After hour thirteen the downwelling begins to trans- 
port the surface water downward, and causes a decrease in the jump, AT, 
strength of a few hundredths of a degree. This combined effect causes 
the deepening of the mixed layer to increase here over that of SI by 
about four meters over the twenty-four hour period. Additionally, it 
decreases the change of temperature of the mixed layer because the layer 
is initially deeper, thereby, reducing mixing, and increasing the dilu- 
tion of the entrained waters, and by causing entrainment of water which 
is not as cool as that of SI, over the same period. These effects are 
apparent in the contours (Figure 16b), and in the T(z) profile (Figure 
16c). The downwelling effect is seen most dramatically in the contours, 
where the downward pulse is evident, along with the described temper- 
ature variations. 

(2) Maximum Downwelling Before Heating (Figure 17 ). Here, 
the downwelling is applied during the last six hours of the mixing. 
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and during the first five hours of the heating cycle. The six hours 
of downwelling during the last six hours of night serve only to deepen 
the layer in the manner previously described. The first five hours of 
downwelling coupled with the heating and wind stress through equation 
(9), produces a noted change in the temperature fields (Figure 17c) 
from that of SI. This coupling leads to an increase of the layer 
temperature at the peak heating period of .008°C for this run over 
that of the previous (ST increase), and over that of the standard, SI. 
After twenty-four hours this difference has decreased to ,005°C as 
compared to the previous run, and by two days the difference in temp- 
eratures is negligible, though the layer depth is slightly shallower 
after two days for this run. By only partially coupling the effects 
of downwelling and heating, a shallower mixed layer has been produced 
with somewhat wanner temperatures. 

The surface temperature evolution and the below layer 
evolution are best revealed by the contoured fields (Figure 17b). It 
is apparent here that another phenomenon is present. As the downwell- 
ing conmences the near layer contours do in fact diverge as expected. 
The deeper contours, however, initially narrow as the pulse forms, 
then widen to the below layer temperature spreading effect described 
previously. This convergence during the peak motion period is a direct 
result of the mathematical treatment of the below layer temperature 
profiles through equation (16). 

(3) Coupled Heating and Downwelling (Figure 18 ). Here 
the downwelling peaks one hour before the maximum surface heat flux. 

The combination of downwelling and wind stress, equation (9) through 
equations (16) and (8), are sufficiently strong to cause the effect of 
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the heat flux, equation (1), to be negligible, and the layer deepens 
to almost three meters more than that of SI by hour eleven. Also, 
the layer temperature has been increased significantly over SI due 
to the downwelling effect on the layer which causes the entrainment 
of less water than previously during the period of maximum surface 
heating. By combining the downwelling with the heating, a maximum 
temperature change of .211°C for SI has been reduced to .162°C for 
this run over the first twenty-four hours. This has also produced a 
mixed layer which is .06°C warmer and three meters deeper than SI 
after only twenty-four hours. This reflects a substantial change in 
the energetics of the system, and may be a significant mechanism for 
the generation of horizontal temperature structure. After the pulse, 
the layer is mixed by the action of the wind stress in the presence 
of a smaller AT— reduced by the downwelling accompanied by a small 
downward surface heat flux. This explains why the deepening is not 
as great after the downwelling pulse as would be expected if one con- 
sidered only the reduced jump strength. 

In combination these three phased runs demonstrate 
the effect of the interaction of the downwelling and the diurnal 
heating cycle. In all cases a relative heating of the mixed layer 
is seen because of the depth increase. The deeper layer reduces the 
effect of the mixing, and reduces the entrainment. As the downwelling 
event couples with the heating, the final mixed layer depth is reduced, 
over that of the uncoupled cases, and greater surface temperatures 
(layer temperatures) are achieved. The rudiments of a pattern in 
surface features which could be induced over a large ocean by simply 
varying the coupling of these two basic effects is now seen. 
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b. Pulsed Downwelling (Case lib) (Figures 19, 20 and 21) 

The initial conditions are those of S2. The magnitudes 
of the downwelling, and the phase relations are the same as the 
previous case. 

(1) 180° Out of Phase with Heating (Figure 19 ). Because 
the downwelling event occurs after the surface heating has stopped (at 
night), the model has gone through the shallowing mode, equations (10) 
and (11), and produced a new surface mixed layer exactly as in S2 dur- 
ing the daylight period. The effect of the downwelling is seen at hour 
fourteen when mixed layer deepening is first observed. The cumulative 
effect of the downwelling is to deepen the layer by about five meters 
more than in S2 in the twenty-four hour period (Figure 19a). The layer 
temperature, due to the downwelling effect, remains ,011°C warmer after 
the period. Destruction of the original layer at 50 m is caused by 
diffusion, and the model's mathematical smoothing in the production of 
new temperature contours due to the downwelling in accordance with 
equation (16). 

(2) Maximum Downwelling Before Heating (Figure 20 ). The 
shallowing of the model through equations (10) and (11) during the 
surface heating period overrides the downwelling, so that by hour five 
this run is essentially the same as S2, with only a slight temperature 
difference. Comparing the contours (Figure 11b and Figure 20b) shows 
that this run is the same as S2 until the downwelling pulse begins at 
hour nineteen, except the early downwelling has deepened the layer 
and moved the cooler water away from the entrainment zone, relative 

to the temperatures in S2. The layer then deepens due to wind stress 
mixing, and downwelling. The net result of this coupled action is a 
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layer about four meters deeper than S2 after twenty-four hours, but 
with a temperature essentially equal to that of S2. 

(3) Coupled Heating and Downwellinq (Figure 21) . Here 
the shallowing mode is still the dominant factor, with the downwelling 
playing only a minor role in that it deepens the newly forming layer, 
and moves the relatively cooler water away from the entrainment zone. 
Because the model is in the shallowing mode, as w e approaches zero, 
the new value of h is dependent on equations (10) and (11), as is the 
temperature of the new mixed layer through equation (1). Since the 
downwelling adjustment is done after the evaluation of the shallowing 
h, the downwelling does deepen the newly formed layer, but no temper- 
ature adjustment is made. This results in a layer that appears to be 
deeper than in S2, and yet maintains the same temperature. The net 
result of these interactions after twenty-four hours is a layer that 
is one meter deeper than in S2, and. .012°C warmer. Here the coupling 
effect was negligible. 

In combination these three runs of the model in the 
downwelling mode with S2 initial conditions showed the opposite pro- 
gression of those runs using the SI initial conditions. The maximum 
variations of temperature and mixed layer depth occurs in the case of 
130° out of phase downwelling with the heating cycle, where a deeper, 
warmer mixed layer is formed during the twenty-four hour period. In 
the case of the SI initial conditions the maximum deviations occurred 
during the coupled run. 

C. Pulsed Upwel ling (Case He) (Figure 22, 23 and 24) 

The initial conditions are those of SI. The phase 
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relationships used are the same as in the pulsed downwelling case, as 
are the magnitudes of the effects. 

(1) 180° Out of Phase with Heating (Figure 22 ). The 
effect here is negligible until hour thirteen, and this run is exactly 
the same as SI until then. At hour thirteen the upwelling begins to 
increase the magnitude of AT. Colder water is moved up toward the 
entrainment zone below the layer, and the layer itself is moved upward, 
as the wind stress acts to deepen the layer by mixing. At hour twenty- 
one it is seen that the combination of the upwelling and the magnitude 
of AT are great enough to actually shallow the layer for that hour. 

The net effect of the upwelling in this case is to 
create a mixed layer that is about 2.6 meters shallower than in SI 
over the twenty-four hour period, and .03°C cooler (Figure 22a). Also, 
there is a maximum 8T in the mixed layer of .241°C, where mixing only 
produced a 8T of ,211°C in SI. This has the potential of producing 
an observable surface temperature effect. 

The effect of contour convergence during the pulse 
mentioned in the downwelling case is seen here again (Figure 22b). 

The overall effect of the upwelling, to converge the contours from 
before to after upwelling, is seen as expected. 

(2) Maximum Upwelling Before Heating (Figure 23 ). In 
this case the upwelling immediately begins to lift the entrainment 
zone against the effect of the wind stress, causing an increase in 
the effect of the mixing, and reducing the water available for 
dilution of the entrained, cold water. In equation (9), w_^ is of 
opposite sign to that of w e , reducing, or changing the sign of dh/dt. 
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This results in cooler layer temperatures over the same period as 
described for SI, and shallower mixed layer depths. Numerically, 
the net result of this effect over the twenty-four hour period is 
to produce a layer that is 2.5 m shallower than that of SI, and ,04°C 
cooler (Figure 23c). The maximum 8T of the layer over the period is 
.240°C, or about the same as in the previous upwelling case. 

(3) Coupled Heating and Upwelling (Figure 24 ). The 
coupling of the upwelling and the heating cycle produce a mixed layer 
that is only 1.9 meters shallower than that of SI in the twenty-four 
hour period, but a layer that is much cooler. The reduction of the 
layer depth during the heating period causes the maximum temperature 
to be .013°C colder, and causes the layer temperature at the end of 
the twenty-four hour period to be .071°C colder than SI, a significant 
change. The maximum 8 T of the layer for this period is ,269°C. 

As in the downwelling case with the shallow initial 
layer conditions, the maximum surface variations in temperature is in 
the coupled upwelling case. 

d. Pulsed Upwelling (Case I Id) (Figure 25 and 26) 

The initial conditions are those of S2. The phase rela- 
tionships used are the same as those in the previous three model runs. 
Magnitudes of the upwelling are the same as in. the previous run. 

( 1 ) 180° Out of Phase with Heating (Figure 25 ) . The 

shallowing proceeds as in S2 for the first six hours, and a new layer 
is formed (Figure 25b). At hour thirteen the upwelling begins and 
several effects are noted. First, the deepening immediately slows, 
relative to that in S2, as the upwelling lifts the layer, and 
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increases the AT at the entrainment zone. Second, the layer temperature 
decreases faster relative to that in S2, due to the shallowing of the 
layer, and entrainment of relatively cooler water into the layer. 

Finally, there is an effect not previously observed. The initial layer, 
which is being diffused in the previous runs is now reinforced during 
the upwelling (Figure 25c), and a two- layer system is formed over the 
first twenty-four hours. 

The net result of this case is a surface layer that is 
2.4 meters shallower and ,022°C cooler than S2 over the period. This 
leads to a slightly increased surface signature over the non-upwelling 
case. Because of the reinforcement of the initial layer there is also 
a potential for a substantially recognizable T(z) signature for this case. 

(2) Maximum Upwelling Before Heating (Figure 26 ). The shal- 
lowing and upwelling combine to make this a relatively uninteresting 
case. Though the large upwelling has some initial effect, the fact 
that w e goes to zero forcing the model to be driven by equations (10) 
and (11) soon overcomes the upwelling, and by hour four this case is 
identical to S2 and the previous run with the exception that the up- 
welling has advected cooler water closer to the entrainment zone. 

The difference in water temperature and mixed layer depth between 
this case and S2 are negligible by hour eighteen. At hour nineteen 
upwelling begins again, increasing the relative effect of the wind 
stress mixing, causing increased entrainment of cool water. The net 
effect of the upwelling and the wind stress mixing over the depth 
involved is to produce a layer that is 2.2 meters shallower than in 
S2, but at only ,0O4°C cooler after twenty-four hours. This would 
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not produce a noticeable effect at the surface. The two-layer type 
system seen in the previous case could, again, provide a signature 
to this sort of upwelling event. 

(3) Coupled Heating and Upwelling (Figure 26) . The 
shallowing prescribed by equations (10) and (11) due to the excessive 
heating tends, initially, to override any effect of the upwelling. 

3y hour six there is only a .4 meter depth difference over that of S2, 
and essentially no temperature difference in the mixed layer. After 
that point, the upwelling begins to combine with the mixing to augment 
the cooling of the layer, and slowing its deepening. 

The net result after twenty-four hours is to produce 
a layer that is about as deep as that of S2, but that is as cool as 
the case of the evening upwelling. This represents a substantial 
energy difference in the runs. Again, a slight surface temperature 
signature could be seen with this run in comparison with the standard. 

3. Case III 

The next test is that of the interaction of mixing with a 
single cycle internal wave. By stepping the phase of an input sinusoid 
through the twenty-four hour progression of the model, a quasi two- 
dimensional motion can be prescribed. 

For a wave of the form 

wave amplitude cC A sin (wt + kx) 
the vertical motion would be of the form 

w(x,z,t) cC rVucos (<ut + kx) (21) 

as seen in Figure 6. By increasing the value of kx, the wave is 
made to travel in the -x direction. 



58 




sanindwv 



59 



Figure 6. Vertical motion-wave form relationship: Vertical motion (dashed line) relative to wave 

(solid line) passing through a given point (x) at time t (hr). Amplitudes are not 
absolute (See equation 21). 



The initial conditions used were those of cases SI and S2. The 
vertical motion was introduced at hour one, and the amplitude was varied 
as t varied in equation (21), to give a wave induced vertical motion. 
This was then run twenty-four times, each time starting at a different 
phase of the wave, to test the interrelation of the wave input with the 
cyclical heat flux. For simplicity, only four major phase relation- 
ships are discussed here, as they will demonstrate the primary effects. 
Also, here, the two-dimensional field will be discussed for the first 
time. 

The wave form of the desired vertical motion does appear in all 
four of the displayed runs. The difference in each run is the phase of 
the vertical motion in relation to the heat flux (daily solar heating 
cycle) . 

a. Shallow Initial Layer (Case Ilia) (Figures 27-33) 

In all four phase relationships presented, maximum downward 
motion at hours 1, 7, 13, and 19, the initial conditions are those of 
SI. The magnitude of the peak upward and downward motions is 36 cm/hr 
at ten feet, as before. The wave is introduced at hour one, and pro- 
gresses through the system through hour twenty-four. 

(1) Peak Upward at Hour 7 (Figure 27 ). In this case the 
initial upward motion, hours one through twelve, the heating and the 
surface wind stress combine to give a slowly deepening layer, which, 
due to the upward advection of the layer and of cooler, below layer 
water, is almost .1°C cooler at its maximum temperature than SI. At 
the end of the upward cycle the layer has deepened about tv/o meters 
less than SI at the same time, and has set up a slightly larger AT at 
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the entrainment zone. As a result, the following downward cycle pro- 
duces a smaller effect in combination with the wind stress in deepen- 
ing than that seen previously. Because of the initial upward motion 
effects on the temperature profile, the combined effect of wind stress 
mixing, surface cooling, and the downward motion produce a layer that, 
after twenty-four hours, is only two meters deeper than SI. Addition- 
ally, the downward effect of the motion on the temperature gradient 
is sufficient only to make up about 60 percent of the deviation in 
temperatures caused by the upwelling motion, as to that of SI. 

The maximum 8T for this case is .245°C. It, also, 
produced the deepest layer of the four runs discussed here, over the 
twenty-four hour period. 

(2) Peak Upward at Hour 1 (Figure 28 ). In this case the 
upward motion was present during the first six hours of the heating 
cycle only. The downward motion cycle t£en progressed for twelve hours, 
followed by six hours of upward motion to complete the twenty-four 
hour sequence. 

During the first six hours the layer deepened to within 
less than one meter of the depth reached by hour six in SI, but the 
layer temperature is seen to be about .01°C cooler. The reversal of 
the motion was too late in the heating cycle to overcome that trend, 
and by hour nine, the time of the maximum surface temperature for SI, 
this experiment produced a surface temperature that is ,013°C cooler. 

The layer depth reached here is within half a meter of the SI depth 
for the same time. The difference is attributed to the initial lift- 
ing of the entrainment zone. 
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As the downward motion continues, along with the sur- 
face cooling, the layer deepens more quickly while the cooling of the 
layer slows as the effect of mixing is reduced. A maximum depth is 
reached two hours after the downward motion has stopped, when the 
upward motion is of sufficient magnitude to overcome the wind stress 
mixing. The layer becomes shallower from its extreme depth of 15.56 
meters to 14.94 meters by hour twenty-four, with associated layer 
cooling due to surface heat loss, and increased entrainment due to 
the increased effect of the wind stress mixing. 

The maximum 8T for this run is ,189°C. The final 
layer depth is within .3 meters of SI for hour twenty-four. The 
variations of depth and temperature during the run are obviously 
much greater than for SI . 

(3) Peak Upward Motion at Hour 19 (Figure 29 ). The fact 
that the downwelling and the heating cycle are in phase, as in Case 

I la, produces, by hour nine, a layer that is two meters deeper and 
,006°C warmer than SI for this same time. The following upward motion 
is sufficient to halt layer deepening at hour sixteen in opposition 
to the mixing, but is not sufficient to cool the layer as much as for 
SI. The final layer temperature is .035°C warmer than that in SI, 
while the final layer depth is about one meter less, though at hour 
sixteen a depth of 14.2 meters was reached. The maximum 8T for tne 
mixed layer in this case was . 181°C. 

(4) Peak Upward Motion at Hour 13 (Figure 30 ). This 
case produced the warmest layer temperature at hour nine, though the 
downward motion was only present during the first six hours of the 
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surface heating cycle. The initial downward motion deepened the layer, 
and thereby reduced mixing in this period, over the SI case, allowing 
the surface layer to be heated to 14.988°C, an increase of ,01°C over 
that in SI. The following upward motion was sufficient, hours eleven 
through fourteen, to hold the mixed layer at a relatively constant 
depth. This allowed a relatively constant, and strong wind stress 
entrainment to combine with surface cooling to rapidly cool the layer. 
The result after twenty-four hours is to produce a layer of about the 
same depth as that of SI, but one that is slightly cooler. 

The 8T here was .237°C, a very noticeable difference 
from SI, though the final product did not vary much in layer temperature 
or layer depth, 

(5) Two Dimensional Data (Figure 31, 32 and 33 ). The two 
dimensional oceanic field was produced as previously discussed. The 
data covered here is that of the evolution of the ocean system at 
hours 6, 12, 18, and 24. 

By hour six the wave has progressed into the oceanic 
field, and its effect below the layer is dramatic, (Figure 31a). It 
has produced a mixed layer that has a mean depth of about 9 meters 
(initial h was 4.5 meters) with a mean temperature of 14.957°C. A 
temperature structure, as could be measured at the surface in x, is 
being formed. The present maximum difference in surface temperatures 
(in x) is ,016°C. 

By hour twelve (Figure 31b) the effect of the wave is 
at its maximum, as seen in the amplitude of the contours below the 
layer. The layer is generally deeper than at hour six with a mean 
depth of about 11.5 meters. The deviation from the mean is much 
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greater for this time, however. The mean of the layer temperature 
has gone to 14.950°C, a slight decrease from hour six, but more 
significant is the deviation of the surface temperatures in the x 
direction. There is now a very definite variation in x of .Q51°C. 

The layer has deepened to a mean depth of about 13 
meters by hour eighteen (Figure 31c), with the effect of the wave 
beginning to be reduced in the system. The mean layer temperature 
has dropped to 14.846°C, a cooling of almost .05°C in the six hours 
between illustrations. The surface temperature deviation has reached 
.091°C. Warm and cold temperature bands are becoming defined in x. 

The single wave has passed through the ocean by hour 
twenty-four, and only the remnant of its passage is seen (Figure 31a). 
In an ocean that is defined by SI conditions everywhere, the twenty- 
four hour product would be an ocean with a mixed layer depth constant 
at 14.65 meters, and mixed layer temperature of 14.767°C. With the 
introduction of a single diurnal wave cycle, a variable condition 
ocean has been produced. The mean depth of this ocean is 15.012 
meters, with local variations of up to +1.73 meters and -1.32 meters, 
with depths at small x being somewhat greater than depths at large x. 
The mean layer temperature is 14.764°C, but more importantly, the 
surface temperature variation is now .095°C, in a clearly defined 
pattern. The pattern that exists is that of a narrow zone of rela- 
tively cool water centered near kx = 5, with a wider zone of relatively 
warmer water centered near kx = 18. Below the layer the isotherms are 
found to be parallel after the passage of this single wave cycle, 
b. Deep Initial Layer (Case 1 1 1 b ) (Figures 34-40) 

The initial conditions for these runs were the same as 
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for S2. The wave cycle passed through the system as in the previous 
runs (Ilia). The wave form is harder to discern in these runs, but 
looking at the contours below the original layer (Figures 34-37) the 
same forms as the previous run are apparent. 

( 1 ) Peak Upward Motion at Hour 7 (Figure 34) . As i n the 
previous deep initial layer cases, the shallowing due to the surface 
heating drives the model in the first six hours. The coincident up- 
ward motion during this time combines with the shallowing and mixing 
to produce a layer that is about 1.5 meters shallower, but at the 
same temperature as in S2, by hour ten. Two noticeable effects during 
this period of upward motion are (i) a slight shallowing of the newly 
formed surface layer, and (ii) the upward transport of cooler below 
layer water, setting up a larger AT at the entrainment zone for this 
case over that of the standard. The downward motion then contributes 
to the deepening of the layer as does the mixing which must initially 
act over a relatively shallow depth. This entrains cool water, form- 
ing a quickly deepening and cooling surface mixed layer. As the down- 
ward motion decreases the cooling slows, and the increase in the depth 
of the layer is seen to be caused mostly by the vertical motion. Of 
all the phases examined this produced the deepest layer (30.34m) and 
has the second coolest layer at the end of the period (15.059°C) in 
relation to the other three cases discussed. 

The 8T for this case is ,115°C. This is about half 
that for the shallow layer case discussed previously. The original 
deep interface can be followed on the contours (Figure 34b) over the 
twenty-four hour period. During the downward motion it, again, is 
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diffusing, while during the upward case there is a slight reinforce- 
ment of the interface. 

(2) Peak Upward Motion at Hour 1 (Figure 35 ). Again it 
is seen that the shallowing mode of the mixed layer model, equations 
(10) and (11), controls the layer development to hour ten. The layer 
depth reached through the shallowing, upward, and downward motion 
combination by hour ten has reached a depth that is about .5 meters 
deeper than that in S2 for the same time. The temperature is the same 
as S2. Also, the downward motion makes entrainment of water into the 
mixed layer due to mixing more difficult, such that cooling after hour 
ten progresses much more slowly than in that of the previous case, and 
that of S2. The net result in the layer depth of this run is a layer 
that is about as deep as that of S2 after the twenty-four hour period. 
The model did not deepen more than S2 due to the dominant upward motion 
during the last six hours. The mixed layer temperature was slightly 
warmer than that of S2 at the end of the run (,071°C). Also, it can 
been seen that the below layer water is somewhat warmer than in S2 for 
the end of the period, due to the combined effects of the heating and 
downward vertical motion timing. 

The 8T for this case is ,092°C. This is less signi- 
ficant than in the previous run. 

(3) Peak Upward Motion at Hour 19 (Figure 36 ). Again, 
the shallowing mode of the model controls the initial ten hour layer 
temperature adjustment. The downward vertical motion, however, drives 
the layer two meters deeper during this time causing reduced mixing 
due to the surface wind stress. The upward motion keeps the layer 
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from deepening as quickly as before after hour twelve. Because of 
the below layer temperature adjustment by the downward motion and 
the initial deepening of the layer, the layer, also, does not cool 
as quickly. This case produced the shallowest final layer depth, 
and the second warmest layer at the end of the period, relative to 
the other three runs. 

The 8T for this case is .103°C, which is about the 
same as the first case which was 180° out of phase with this run. 

(4) Peak Upward Motion at Hour 13 (Figure 37 ). The 
coincident shallowing, wind stress mixing, heat flux, and cycled 
vertical motion over the first ten hours produce a surface mixed 
layer that has the same temperature as S2, and almost the same layer 
depth (to within less than .5 meters). Below the layer, however, 
the upward motion has made relatively cooler water available by hour 
ten (in comparison with S2), and the layer, though not rapidly deep- 
ening due to the continued upward motion, cools quickly as large 
quantities of water are entrained over a relatively shallow depth, 
a trend which continues to the end of the period. The net result is 
a mixed layer that is 1.5 meters deeper than that of S2, but is also 
the coolest layer of the four runs. Because of the below layer temp- 
erature profiles which are set up by the upward motion for the mid- 
twelve hours, the continued run of this model over a second day gives 
continued rapid cooling of the layer due to increased entrainment of 
cooler, upward advected water, and to increased effects of mixing. 

A noticeable deepening does take place at the end of the period due 
to the downward motion, which accounts for the fact that this run 
ended deeper than S2. 
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The 8T for this run is . 125°C. This is the largest 
of the four runs, but it is still only half of the maximum 8T seen for 
the shallow initial layer case. 

(5) Two-Dimensional Data (Figure 38, 39 and 40) . The two- 
dimensional oceanic field is produced, here, as previously discussed, 
for the single wave cycle. The data covered is for hours 6, 12, 18, 
and 24. A common feature is the varying below layer temperature struc- 
ture set up because of the introduction of the initial deep layer. 
Because of this, the changes in the thermocline below the initial deep 
interface best show the effect of the wave progression through the 
system. 

As in the shallow case, by hour six, the wave has 
progressed into the oceanic system producing some horizontal sea sur- 
face temperature structure. However, the effect on the surface 
temperature is almost negligible, in that the variability of the 
surface temperature in x is only ,001°C. This is due to the shallow- 
ing mode response of the model over the first six hours. Even with 
this, however, the depth structure has varied (Figure 38a) so that 
there has been a change from the constant structure that would be 
seen in S2 of 11.767 meters for the layer depth, to a layer of mean 
depth of 11.769 meters t .35 meters, indicating that an adjustment 
due to the wave is taking place. 

By hour twelve much more temperature and layer ( h(x) ) 
variation is beginning to take shape. At hour twelve the case without 
the vertical motion has a layer depth of 15.037 meters, and a layer 
temperature of 15.159°C everywhere. In this case, where the maximum 
effect of the wave has been reached, the layer depth has a mean of 
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15.231 meters with a deviation of t 2.5 meters. The variability in 
the surface temperature in x has reached .01 °C. Though the variabil- 
ity in the surface temperature is not very large, a pattern is again 
beginning to emerge in the x direction. At small x cooler water is 
forming, while at large x a band of warmer water is emerging. 

The layer has continued to deepen to hour eighteen. 
The layer depth has increased to a mean of 19.46 m, t 3.4 m, and the 
variation in the surface temperature structure in x has reached ,03°C. 
The warm-cold pattern seen at hour six is repeated and strengthened. 
Though the temperature difference between warm and cool water is not 
as great as in the shallow case, the banding of the patches is much 
more distinct here (Figure 33, Figure 31). 

At the end of the period the wave has, again, been 
removed from the system, and only the remnant of its passage is seen. 
The final layer depth is deeper than in the previous hours, and much 
deeper than in the shallow case. The variability of h in x is also 
much greater here than in the shallow case. The surface temperature 
patchiness is well-defined, but the variability of the temperature in 
x is only ,035°C. 

Comparing the two cases, shallow and deep initial 
layers, it is seen that for the same amplitude wave the model produces 
two distinct ocean types. For the shallow case an ocean is formed 
that has a relatively shallow layer with a large temperature varia- 
tion in the x direction, but with that temperature variation primarily 
in wide bands. For the deep case a deep surface layer is formed with 
smaller variability in the surface temperature. However, the temp-, 
erature structure in x is much better defined. Also in the deep case. 
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the patches of warm and cold water are much closer to each other (in x). 
This could suggest a mechanism for generating or sharpening SST fronts 
through surface heating, and internal solitary wave passage. 

4. Case IV 

It was noted in Case III that after twenty-four hours the wave 
form had passed through the ocean system, and had left behind only the 
residue of its effect on the system, over and above the standard effects. 
This can be compared to a solitary wave passing through an ocean system. 
While this has theoretical advantages, it is not a very realistic prob- 
lem, in that wave forms seldom take on this sort of motion. One pos- 
sible exception may be soli tons generated at extreme tidal stages 
(Apel , 1982). 

The next progression, therefore is to look at a continuous wave 
system that passes through the ocean with constant form, until the run 
is terminated. This is done by simply modifying the mathematics of 
the previous case to include initialization of the temperature and 
salinity profiles with the wave form according to the following. 

We know that 



w(x,z,t) = AAz_ cos(<ut + kx) 



and 
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therefore 



For initial displacement 

8h(x,z t O) = 8h 0 

therefore 

8h 0 = -AAz sin(kx) 

wD (22) 

Applying equation (22) to form the initialized layer depth, and com- 
bining this with equation (16) to initialize the temperature profile, 
a continuous wave field can now be passed through the model ocean 
(Figure 7). 

From this point on the discussion will be confined to the deep 
initial layer case, $2 standard, and variations thereof. This is done 
because all the primary effects discussed to this point appear in both 
the shallow and deep cases, as do the effects to be discussed in vary- 
ing degrees, and it is not necessary to continuously reiterate this 
discussion. 

a. Diurnal Wave (Case IVa) (Figures 41, 42 and 43) 

Because of the overriding effect of the shallowing mode of 
the model in the first six hours of the runs, the results produced 
here were exactly the same as those discussed in Case 1 1 1 b , for the 
phase relationships. This was of course expected from the theory. 

(1) Two-Dimensional Data (Figure 44, 45 and 46 ). Again, 
Decause of the dominance of the shallowing mode, the hour 6, 12, 18, 
and 24 two-dimensional fields produced here are the same as those of 
Case 1 1 1 b for the area of the entrainment zone, and the mixed layer. 
3elow the layer, however, the wave form is now seen to progress through 
the system in a continuous mode. This is best seen in the contours 
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Figure 7. Diurnal wave relationship (of continuous diurnal waveform to model location (x). 
Amplitude variable.) 



in which the amplitude of the system remains a constant, but the wave 
form is shifted with time. 

V 

It can be seen that the patch of warmer waters forms 
(Figure 44) and takes a phase that barely leads the peak of the wave 
system. The entrainment zone ( h(x) ) now takes the form of the 
passing wave system, as would be expected, and yet was not obvious in 
the previous case. The divergence and convergence of below layer 
fields is also obvious here, where previously it was not particularly 
obvious. 

b. Semi-diurnal Wave (Case IVb) (Figures 47-50) 

A dominant internal wave system of the temperate oceans 
is the semi-diurnal tidal cycle. South of 30° N there is also a 
strong diurnal tidal cycle. A semi-diurnal tidal field can be pro- 
duced for testing in this model by simply adjusting the frequency of 
the prescribed w for equation (21), and setting the initial profiles 
with the appropriate frequency of equation (22) and (17). 

The amplitude of the motion remains the same as the 
diurnal cycle case. Secause the cycle of wave mixing interactions 
repeats itself after a 2 ir phase shift, which occurs after 12 hours 
in this case, only the cases of six and twelve hour phase shifts 
will be discussed. The 18 and 24 hour phase shifts are merely repeats. 

(1) Six Hour Shift (Figure 47a and b ). Figure 8 shows 

4 

the relationship of the heating cycle to the wave motion in this case. 
The vertical motion is seen to immediately combine with the shallowing 
mode of the model through the initialization to give a somewhat 
shallower initial layer. The shallowing mode of the model holds the 
layer temperature the same as in S2 to hour six, while the vertical 
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Figure 8. Phase relationship of Case IVb (six hour shift) of semi-diurnal waveform to diurnal 
surface heating cycle over time (hrs). Amplitudes are not proportional. 



motion varies the layer depth, and alters the below layer temperature 
structure. By hour six a half-cycle of the wave has passed through 
the model point, and a layer has been produced that is almost one 
meter deeper than S2 at this time, but which has the same layer temp- 
erature. Continued but decreasing downward motion combines with the 
downward heat flux until hour nine to slowly deepen the layer, but 
the layer continues to warm. This process continues through hour ten, 
though upward motion has begun because the heating and entrainment of 
relatively warm water are greater than the effect of the resistance to 
deepening due to the surface heat flux, and the small upward motion. 

By hour eleven the upward motion dominates, and the layer shallows 
while it cools slightly due to the increased entrainment. At hour 
twelve the layer is as shallow as it will ever get, and the first wave 
has passed. Now decreased upward motion followed by downward motion, 
combined with surface cooling and mixing through a progressively cooler 
surface mixed layer quickly increases the layer depth to a maximum of 
26.7 meters at hour twenty-one. This is deeper than S2 gets in the 
entire twenty-four hour period. Yet, the layer is not as cool as in 
the S2 case. By hour twenty-one the upward motion is again dominant 
and the layer shallows. 

(2) Twelve Hour Shift (Figure 47c and d ). Figure 9 shows 
the relationship of the heating cycle with the wave in this case. The 
initialization combines with the vertical motion and shallowing mode 
to produce an initial layer that is deeper than that in S2, but soon 
becomes shallower, as the new surface layer is established. By hour 
six a layer of the same temperature as S2, but .72 meters shallower 
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Figure 9. Phase relationship of Case IVb (twelve hour shift). 



has been formed. The net result of the cyclical interactions in this 
case is a layer that is 1.5 meters deeper than in S2, but which is at 
about the same temperature. The difference in the layer temperature 
between this run, and the previous is .004°C. However, the layer 
depth is 2.2 m different. 

(3) Two-Dimensional Data (Figure 48, 49 and 50 ). At hour 
six a very discernable depth structure has been set up at the bottom 
of the mixed layer which reflects the structure of the passing wave. 

Above the entrainment zone the temperature structure, as produced by 
the shallowing mode, is homogeneous. Below the layer the structure is 
highly varied in x reflecting the combined effects of the vertical 
motion to this time. The mean layer depth is 11.75 meters, while that 
of S2 is 11.77 m, showing the domination of the shallowing to this point. 
A variation of h in x of ± .75 meters is already exhibited. 

By hour twelve the mixed layer depth structure in x, 
and the layer temperature variability in x have evolved considerably 
over that of hour six. At the entrainment zone the structure is approx- 
imately sinusoidal, as expected. Over the last six hours a large 
amount of water from below the layer has been entrained while surface 
heating has occurred, and the mixed layer has heated and cooled 
accordingly to begin to set up the pattern of warm and cold bands of 
temperature in x. The mean depth of the layer is now 15.21 meters 
t 2.50 meters, while that for S2 is 15.04 meters. The temperature 
variation in x as seen at the surface is now ,012°C. 

The structure of temperature and layer depth profile 
has been further accented by hour eighteen. The layer has continued 
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to deepen, and the fluctuation in x increased. The banding of the 
in layer temperature pattern in x has become more defined. 

After twenty- four hours a layer has been produced 
that has a depth of 26.14 m, t 3.3 m. The surface temperature pattern 
in x is quite clearly defined, though the variation in temperature is 
only .015°C. 

5. Case V (Superimposed Waves) (Figure 51, 52 and 53 ) 

Since oceanic internal waves are almost never simple sinusoids, 
one more progression is considered. Spectral analysis of the MILE 
data (Figure 10) at station P (55° N) shows the dominant wave field to 
be the semi-diurnal tide superimposed by, at least, the 6-7 hour har- 
monic of the tidal wave. Internal wave theory indicates that the 
period of the internal waves must be less than or equal to the inertial 
period. North of approximately 30° N, no diurnal internal waves may 
propagate, and therefore, they are not seen in the MILE data. South 
of this, however, diurnal and semi-diurnal tidal period waves may 
propagate. Therefore, since previous discussions have been concerned 
with diurnal and semi-diurnal waves, and there is justification for 
their presence, these waves will be superimposed to form a slightly 
more complicated wave system. 

The combined wave is directed through the two-dimensional 
ocean which has the initial conditions of S2, from right to left 
(large kx to small kx). The amplitude of the motion is such that 
both waves peak at the same amplitude at their peak motions, though 
one has half the period of the other. Both waves start with an 
amplitude of vertical motion in phase at time one. 
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Figure 10, MILE spectral analysis (of temperature variation at 50 in (□) , 47 m (+), and 32 in (C$). 



The superimposed waves reflect all of the interactions discussed 
in the previous sections. Further, this simple addition of two wave- 
forms has developed several new features. The new features can be 
traced in their development from hour one through hour twenty- four, 
where they become most evident. Since this is the case, the discussion 
will be limited to the case of the two-dimensional ocean at hour twenty- 
four. 

By this time a complex structure has emerged for the two-dimen- 
sional ocean, both in temperature and in profile of the layer depth. 

An ocean has been formed that has a mean layer depth of 26.78 m, but 
this depth has a variation that is not symmetrical from small to large x. 
At small x the variation from the mean depth is seen to be a maximum of 
-5.45 m, whi\e at large x the layer depth varies to a depth deviation 
from the mean of +7.83 m. An ocean has been formed that has a deep 
layer at one end, and a shallow layer at the other. Further, a devia- 
tion of the surface temperature in the x direction of .042°C is seen. 

The pattern of the warm and cold bands again is seen, but this also is 
no longer symmetrical. The temperature pattern, as seen at the surface, 
is a cold band beginning at small kx, progressing through a weak gradi- 
ent to a warm band centered near kx = 15. A strong gradient is then 
encountered leading to a cold band centered about kx = 21. Through the 
introduction of a simple two wave system, an extremely diverse ocean 
has been produced. 

Shifting the phase relationships of the two internal waves, and 
running the model again, resulted in similar patterns. Positions of 
the temperature maxima and minima in the layer, and the position of 
the depth variations shift, but their relative positions are maintained. 
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The temperature variation in the x direction as seen at the surface 
does not change with the phase changes. 
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IV. CONCLUSIONS AND RECOMMENDATIONS 



The analysis of the mixing response to the prescribed vertical 
motion shows that vertical advection can be successfully introduced 
into this one-dimensional mixed layer model without numerical diffi- 
culties. This is important in that it indicates that a single station 
forecasting model can be altered to include the variability of long 
internal wave motion. This step is the first of a progressive devel- 
opment of a single station forecasting model that may in the future 
include other local advective effects. 

Using an ensemble of one-dimensional stations, a two-dimensional 
mixing response to prescribed internal wave motion was simulated. 
Although the initial temperature field (without wave) and the surface 
boundary conditions were both horizontally homogeneous, a significant 
wave-induced horizontal structure was created in the T(x,z) field, doth 
linear and non-linear effects have been demonstrated in the simple 
single wave and superimposed wave cases. Horizontal variability has 
been shown, both in mixed layer depth and in mixed layer temperature 
structure. The mixed layer temperature develops a corrollary structure 
of alternating warm and cold patches. In the case of the superimposed 
waves of different frequencies, this patchiness pattern can be complex, 
and can give rise to horizontal gradients which could be frontogenetic. 

There is an important relationship between the phase of internal 
motion and the cyclical diurnal heat flux boundary condition. This is 
clearly demonstrated in the numerical experiments in which the phase 
of the input pulse or wave system has been varied. 
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In the case of superimposed diurnal and semi-diurnal waves, the 
horizontal variation in temperature in the mixed layer was about .04°C 
after only a single diurnal cycle. In the same period, the mixed layer 
depth varied by about 13 m. The wave dispersion relation for a two layer 
ocean at 25° N gives a horizontal wave length scale of from 10 to 100 Km, 
depending upon stratification. This wave-mixing interaction may, there- 
fore, be an important mechanism for the generation of sea surface 
patchiness as found by Butler (1981). 

Further study must be done in this area. Clearly, more complicated 
waveforms will produce more complex responses. The model should be 
tested with salinity changes, including evaporation and precipitation. 
Tests should be conducted on data sets such as those of MILE. Compar- 
ison of the mixed layer model results using MILE boundary conditions, 
with and without vertical advection, as observed in spectrum analysis, 
should further demonstrate the importance of addition of vertical 
motion to mixed layer models. 

Finally, the wave systems should be tested in an interactive mode 
with the mixing of the model. Energy exchanges, as prescribed by wave 
theory and the studies of Linden (1975), and Kantha et al. (1977) 
should be allowed, and the relative significance of the different 
mechanism should be evaluated. Internal vertical motion is an impor- 
tant dynamic variable, and should be treated as such in the future of 
mixed layer modelling. 
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Figure 17. Same as figure (1), except 
for Case Ha with peak at hour 24. 



13.00 



: 3.50 

~ E ” P . 



14.00 14.50 15.00 

'RELATIVE FIRST) 



15.50 16.00 16.50 



90 




91 



OAT HR H 

11 49.5 

1 2 18.1 
1 3 16.2 

1 4 13.0 

1 5 12.1 

1 6 11.8 

1 7 12.1 

1 8 12.5 

L 9 13.0 

1 10 13.5 

1 11 14.2 

1 12 15.0 

1 13 16.0 

1 14 17.2 

1 15 18.6 

1 16 20.0 
1 17 21.5 

1 18 23.1 

1 19 24.7 

1 20 26.2 
1 21 27.7 

1 22 23.9 

1 23 30.0 

1 24 30.3 



T S 

15.001 34.000 
15.009 34.000 
15.024 34.000 
15.047 34.000 
15.075 34.000 
15.105 34.000 
15.132 34.000 
15.154 34. OOC 
15.168 34.000 
15.175 34.000 
15.172 34.000 
15.160 34.000 
15.149 34.000 
15.139 34.000 
15.130 34. COO 
15.122 34.000 
15.114 34.000 
15.107 34.000 
15.101 34.000 
15.096 34.000 
15.090 34.000 
15.086 34.000 
15.081 34. OOC 
15.077 34.000 



TAU QNET 

0.50 - 55.7 
0.50 -201.5 
0.50 -326.8 
0.50 -423.0 
0.50 -483.5 
0.50 -504.2 
0.50 -483.7 
0.50 -423.5 
0.50 -327.5 
0.50 -202.4 
0.50 - 56.6 
0.50 99.9 
0.50 100.8 
0.50 100.8 
0.50 100.8 
0.50 100.8 
0.50 100.8 
0.50 100.8 
0.50 100.8 
0.50 100.8 
0.50 100.8 
C.50 100.8 
0.50 100.8 
0.50 100.8 



WOO OH 

0.0 O.C 

0.0 0.0 

0.0 0.0 

0.0 0.0 

0.0 0.0 

0.0 0.0 

0.0 0.0 

0.0 0.0 

0.0 0.0 

0.0 0.0 

0.0 0.0 

0.0 0.0 

0.00932 14.869 

0.01800 30.739 

0.02546 46.756 

0.03118 61.334 

0.03477 73.399 

0.03600 81.528 

0.03477 84.247 

0.03118 80.442 

0.02546 69.511 

0.01800 51.613 

O.OC932 27.799 

0.0 0.0 



2 1 31.1 

2 2 2 1.0 

2 3 15.4 

2 4 13.1 

2 5 12.0 

2 6 11.7 

2 - 12.0 
2 9 12.4 

2 9 12.9 

2 10 13.4 

2 11 14.1 

2 12 14.9 

2 13 15.9 

2 14 17.0 

2 15 19. 3 

2 16 19.8 

2 17 21.3 

2 IS 22.9 
2 19 24.6 

2 20 26.2 
2 21 27.6 

2 22 28.9 

2 23 30.0 

2 24 30.8 



15.077 34.000 
15.085 34.000 
15.101 34.000 
15.127 34.000 

15.152 34. COO 
15.181 34.000 
15.208 34.000 
15.231 34.000 
15.245 34.000 
15.251 34.0CG 
15.248 34.000 
15.236 34.000 
15.225 34.000 
15.215 34.000 
15.206 34. COO 
15.198 34.000 
15.190 34.000 
15.183 34.000 
15.177 34.000 
15.172 34.000 
15.166 34.000 
15.162 34.000 
15.157 34.COO 

15.153 34.000 



0.50 -55.7 
0.50 -201.5 
0.50 -326.8 
C.50 -423.0 
0.50 -483.5 
0.50 -504.2 
0.50 -483.7 
0.50 -423.5 
0.50 -327.5 
0.50 -202.4 
0.50 - 56.6 
0.50 99.9 
0.50 100.8 
C.50 IOC. 8 
0.50 100.8 
0.50 100.6 
0.50 100.8 
0.50 100.8 
0.50 ICO. 8 
C.50 100.9 
0.50 100.8 
0.50 100.8 
0.50 ICC. 8 
C.50 100.8 



0. 0 
0.0 
0.0 
0.0 
C. 0 
0.0 
0.0 
O.C 
0.0 
0.0 
0.0 
0.0 

0.00932 

0.01800 

0.02546 

0.03118 

0.03477 

0.03600 

0.03477 

0.03118 

0.02546 

0.01800 

0.00932 

0.0 



0.0 
0.0 
0.0 
0. 0 
0.0 



0. 0 



14. 752 
30.358 
46.032 
60.690 
72 . 794 
81.054 
83.923 
80 .282 
69.477 
51.644 
27.833 
0 . 0 



(a) 

3 

O 

nl 




§! 

SH ! 



5.3000 
,5.7090 
15.0240 
15.0470 
1S.Q7S0 
15. 1040 

is.mo 

5. 1 5 90 
15. 690 
5.1740 
5.1710 
*5. 1590 
15. ! «80 
: 5 . : 380 
5.1290 
15.1210 
5. 11<*0 
.5.1370 
5.1010 
5.C9S0 
.5. 0930 
.5.3950 
15.0910 
•5.0760 



-49.5152 
19. 1073 
16.1 926 
3.0190 
12.0605 
11.7670 
12. 1436 
.2.4381 
12.9925 
13.5208 
14.2001 
15.017? 
16.0323 
17.2313 
18.5022 
19.9913 
21. -770 
23.357! 
24.5510 
25.2063 
27.5SS9 
29.3314 
29.9735 
30.7566 



(c) 



Figure 19. 5a«e as figure Cl), except 
tor Case I lb with peak at hour 18. 
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15.154 
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0.0 


0.0 


l 


13 
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15.149 
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14 
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IS 
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16 
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15.121 


34 . 000 


0.50 
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l 


17 


19.3 


15.113 
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0.50 


100.8 
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i 


18 
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15.105 
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O.SO 
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0.0 


0.0 


l 


19 


21.3 


15.098 


34.000 
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0.00932 


19.714 


l 


20 


22.5 


IS. 09 1 


34.000 
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100.8 


0.01800 


40.192 


l 


21 


24.0 


15.385 


34.000 
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loc. a 


C. 02546 


60. 364 


1 


22 


25 . 7 


15.079 


34.000 


O.SO 


100.8 


0.03118 


78.928 


l 


23 


27.6 


15.074 


34.000 


0.50 


100. 8 


0.03477 


94.207 


l 


24 


29.5 


15.069 


34.000 


O.SO 


100.8 


0.03600 


104.386 


j 


1 


31.0 


15.070 


34.000 


0.50 


- 55.7 


0.03477 


106.008 


2 


2 


2 X . 7 


15.077 


34.COO 


0.50 


-201.5 


0.03118 


66.710 


2 


3 


15.6 


15.093 


34.300 
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-326.8 


0.02546 


39.216 


2 


4 


13.3 
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34.000 


0.50 


-423.0 


0.01800 


23.782 


2 


5 


12.1 


15.144 


34.000 


0.50 


-463.5 


0.009)2 


11.182 


2 


6 
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15.173 


34.000 


0.50 


-504.2 


0.0 


O.C 


2 


7 


12.0 


15.201 


34.000 


0.5C 


-433 . 7 


0.0 


0.0 


2 


8 


12. 4 


IS. 22 3 


34.000 


0.50 


-423.5 


0.0 


0.0 


2 


9 
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15 .236 


34.000 


0.50 


-327.5 


0.0 


0.0 


2 


13 


13.5 


15.244 


34.000 


0.50 


-202.4 


0.0 


0.0 


2 


IX 


14.2 


15.2 41 


34.000 


0.50 


- 56.6 


0.0 


0.0 
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75.0 


15.229 


34. COO 


C.SO 


99.9 


0.0 


c.o 


2 


13 


15.8 


15. 218 


34.000 


0.50 


100.8 


0.0 


0.0 


2 


14 


16.7 


IS. 208 


34.000 
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0.0 


0.0 


2 


15 


17.5 


15.198 


34.CC0 


0.50 


100.8 


0.0 
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2 


16 


18.3 


15.139 


34.000 


0.50 


ico. a 


0.0 


0.0 


2 


17 


19.2 


15.181 


34.300 


0.50 


ioo. a 


c.o 


0.0 
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18 


20.1 


15.173 


34.000 


0.50 


100.8 


0.0 


0.0 


2 


19 


21.2 


IS. 166 


34 .COO 


0.50 


ICC .8 


C .00932 


19.677 


2 


20 


22 . 5 


15.159 


34.000 


C.SO 


ico. a 


0.01800 


40.216 


2 


21 


24.1 


15. 153 


34.000 


O.SO 


100.8 


0.02546 


60.511 


; 


22 


2 S . 3 


15.147 


34.000 


0.50 


100.8 


0.03118 


79.246 


2 


2) 


27. 7 


15.142 


34.000 


0.50 


100.8 


0.03477 


94.597 


: 


24 


2 9.6 


25.13- 


34.000 


0.50 


100.8 


0.0360C 


104.837 










(a) 













TIME (HR) 



_arc« ’E-p. OEPTh 



JS.OCOO 
15.0090 
1 5. 02*0 
:5. 34^0 
15.0750 

is. :o«»o 

4 15.1)10 

15. 15-0 
15. 1600 
15.1 NO 
15.H20 
15. 1500 
IS. U50 
15. 1 390 
i 15.1290 

15.1210 
15. 1120 
15. 1050 
15.0970 
15.0910 
:5.09«*0 
j 15.0793 



51.267) 
18. '•99! 
16. -.907 
15.1960 
12.1)30 
11.7510 
12. 1393 
12. *952 
13.3005 
13.5514 

:4.2h:9 
15.0344 
15.3397 
16.7917 
17.6410 
9.-385 
19.3379 
20.1931 
21.2566 
22.5307 
24.0162 
25. 7130 
27.5654 
29.5215 



(C) 



| 



Figure 20. Sa»e as figure (1), except 
for Case lib with peak at hour 24. 




12.00 12.50 



13.00 



13.50 u700 14.50 15.00 
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16. CO 16.50 
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0.03477 


45. 7S0 


1 


5 


12.4 


IS .076 


34.000 


0.50 
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-504.2 


0.03477 
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12.6 


15.133 
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0.50 


-483,7 


0.03118 


38. 787 
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9 


13.3 
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34.000 
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-423.5 
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-327.5 
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19 
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15.107 
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0.0 


l 


20 
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38.240 
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9 
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IS. 232 


34.000 


o.so 


-423.5 


0.02S46 


33.120 


2 


9 


14.0 


15.247 


34.000 


O.SC 


-327. S 


0.0180C 


24.907 


2 


10 


14.7 


1S.2S3 


34.000 


0.50 


-202.4 


0.0C932 


13.616 


2 


11 


IS. 4 


1S.2S1 


34.000 


o.so 


- 56.6 


0.0 


0.0 


2 


12 


36.2 


IS. 241 


34.000 


o.so 


99. 9 


0.0 


0.0 


2 


13 


17.1 


IS. 231 
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0 . 50 


100.8 


0.0 


0.0 
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14 


17 . 9 


IS. 222 


34.000 


o.so 


100.8 


0.0 
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15 
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IS. 213 
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100.8 


0.0 


0.0 
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16 
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34.000 


o.so 


100.8 


0.0 
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17 


20.3 


IS. 197 
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o.so 


100.8 


0.0 


0.0 
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18 


21.1 
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34.000 


o.so 


100.8 


0.0 
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22.0 
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34.000 


o.so 


100.8 


0.0 


0.0 


2 


20 


22.8 
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o.so 


100.8 


0.0 


0.0 


2 


21 


23. 7 
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34.000 


o.so 


100.8 


0.0 


0.0 


2 


22 


24.6 


IS. 16S 


34.000 


o.so 


100.8 


0.0 


0.0 


2 


23 


25. S 


IS. 159 


34.000 


o.so 


100.8 


0.0 


0.0 


2 


24 


26.7 


IS.1S4 


34.000 


o.so 


100.8 


0.009 32 


24.791 



(a) 



.5.0003 
15. 0093 
:s.C?i3 
5.3473 
5.0750 
5. 3S0 
.5.1 320 
.5.1540 
•5. .590 
15.1750 
15.1 7 JO 
15.1530 
5. .533 
5. 1440 
5. 1 3S0 
.5. .2*3 
.5. 1 20C 
IS. .30 
5.1 060 
•5. 1C00 
s.:9«o 
5. 0990 
5.0920 
5.0770 



S0.41HS 
9.4797 
15.5906 
1 3. 3«6* 
.2.4227 
12.1067 
2.5361 
1 3 3217 
.4.0359 
4.9357 
S.6259 
15.4615 
7.2'91 
5.0945 
.9. 3090 
19. 7240 
23.5437 
21 . J563 
22. 1325 
23.3234 
23.5532 
24.6944 
25.5522 
26.6794 



(c) 



Figure 21. S ue as figure (1), except 
for Case lib with peak at hour $. 





TIME fHR) 

(b) 



. 2 .:a j.so : 3 .oo n.so m.oo w.sa is.oo is. so 16.00 t 6 .so 

TE“P. 'RELATIVE TO FIRST) 
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OtPIH IMJ 

60.00 bJ.OO 46. OU 39.00 -32.00 25.00 1 6^0 II . UO 4.00 3 y 00 




12.00 :2.S0 13.00 13.50 W.CO U.SO 15.00 IS. 50 16.00 16.50 

' r E M P . :PELfl T IVE *0 FIRST! 
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96 



1 



OAT HR H 



T 



S TAO QNST WDD 



OR 



1 1 47.8 

1 2 17.7 

1 3 IS. 9 

1 4 12.9 

1 5 12.0 

l 6 11. 8 

1 7 12.1 

1 8 12.5 

1 9 12.9 

1 10 13.5 

L 11 14.1 

1 12 1 S .0 

1 13 15.8 

1 14 16.6 

1 IS 17.4 
1 16 18.2 
1 17 19.0 

1 19 19.9 

1 19 20.6 

1 20 21.1 
1 21 21. S 

1 22 21.8 

1 23 22.1 

1 24 22. 5 



15.001 34.000 
15.009 34.000 
15.024 34.000 
IS. 047 34.000 

15.075 34.000 
15.104 34.000 
IS. 131 34.000 
15.153 34.000 
15.168 34.000 
15.174 34.000 
15.171 34.000 
15 . IS 9 34.000 
15.148 34.000 
15.138 34.000 
IS. 128 34.000 
IS. 119 34.000 
15.111 34.000 
15.103 34.000 
15.095 34.000 
15.088 34.000 
IS. 082 34.000 

15.075 34.000 
IS. 068 34.000 
IS. 062 34.000 



0.50 - S5.7 
0.50 -201.5 
0.50 -326.8 
0.50 -423.0 
0.50 -483.5 
0.50 -504.2 
0.50 -483.7 
0.50 -423. S 
0.50 -327.5 
0.50 -202.4 
O.SO - 56.6 
0.50 99.9 
0.50 100.3 
O.SO 100.8 
0.50 100.8 
0.50 100.8 
0.50 100.8 
0.50 100.8 
0.50 100.8 
0.50 100.8 
0.50 100.3 
O.SO 100.8 
0.50 100.8 
0.50 100.8 



-.03477 
-.03118 
-.02546 
-.01800 
- . 00932 
-.00000 



0.0 



-.00932 

-.01800 

-.02546 

-.03118 

-.03477 

-.03600 



-169.223 

- 56.143 

- 40.938 

- 23.344 

- 11.218 

- 0.000 

0.0 



0.0 

0.0 



0. 0 



- 19.256 

- 38.331 

- 55.477 

- 69.136 

- 78.296 

- 82.402 



2 1 22.4 

2 2 21.9 

2 3 14.3 

2 4 13.0 

2 5 11.9 

2 6 11.3 

2 7 12.1 

2 9 12.4 

2 9 12.9 

2 10 13.4 

2 11 14.1 

2 12 14.9 

2 13 15.7 

2 14 16.5 

2 15 17.3 

2 16 18.2 
2 17 19.0 

2 19 19.3 

2 19 20.4 

2 20 20.7 

2 21 20.8 
2 22 20.8 
2 23 20.7 

2 24 20.6 



15.062 34.000 
15.069 34.000 
15.085 34.000 
15.108 34.000 
15.136 34.000 
IS. 165 34.000 
15.192 34.000 
IS. 214 34.000 
15.229 34.000 
IS. 235 34.000 
15.232 34.000 
15.220 34.000 
15.209 34.000 
15.199 34.000 
15.189 34.000 
15.180 34.000 
15.172 34.000 
15.163 34.000 
15.156 34.000 
15.149 34.000 
15.142 34.000 
15.134 34.000 
15.127 34.000 
15.120 34.000 



O.SO - 5S.7 
0.50 -201.5 
0.50 -326.8 
0.50 -423.0 
O.SO -493.5 
O.SO -504.2 
0.50 -483.7 
0.50 -423.7 
O.SO -327.5 
O.SO -202.4 
0.50 - 56.6 
0.50 99.9 
0.50 100.8 
0.50 100.8 
0.50 100.8 
O.SO 100.8 
O.SO 100.8 
0.50 100.8 
0.50 100.8 
0.50 100.8 
O.SO 100.8 
O.SO 100.8 
O.SO 100.8 
0.50 100.8 



-.03477 
-.03118 
-.02546 
-. 01800 
-.00932 
-.00000 



-.00932 
-.01800 
- .02546 
-.03118 
-.03477 
-.03600 



79.400 
69.231 
38.163 
23.687 
11. 136 
0.000 



0.0 

0.0 



0.0 



0.0 



0.0 

19.06S 

37.614 

53.606 

63.729 

73.101 

75.372 



(a) 

o 

o 





97 







98 



1 



OAT 


HR 


8 


T 


S 


TAO 


QNET 


W DO 


oa 


1 


1 


6 . 3 


14.956 


34.000 


0.50 


- 3.5 


0.0 


0.0 


1 


2 


7.3 


14.939 


34.000 


0.50 


- 100.8 


0.0 


0.0 


1 


3 


9.0 


14.935 


34.000 


0.50 


- 184.3 


0.0 


0.0 


1 


4 


a. s 


14.939 


34.000 


0.50 


- 248.4 


0.0 


0.0 


1 


5 


9.9 


14.949 


34 .000 


0.50 


- 288.7 


0.0 


0.0 


1 


6 


9 . 2 


14.958 


34.000 


0.50 


- 302.5 


0.0 


0.0 


1 


7 


9.5 


14.969 


34.000 


0.50 


- 288.9 


0.0 


0.0 


1 


8 


9.9 


14.976 


34.000 


0 . 50 


- 248 . 7 


0.0 


0.0 


1 


9 


10. 1 


14.978 


34.000 


0.50 


- 184 . 7 


0.0 


0.0 


1 


10 


10.4 


14 .974 


34.000 


0.50 


- 101.3 


0.0 


0.0 


1 


11 


10.7 


14. 963 


34.000 


0. 50 « 


- 4.1 


0.0 


0.0 


1 


12 


11.1 


14 .943 


34.000 


0.50 


100.2 


0.0 


0.0 


1 


13 


11.3 


14.926 


34.000 


0.50 


100.8 


-.00932 


- 10.604 


1 


14 


11.5 


14.909 


34.000 


0.50 


100.8 


-.01800 


- 20.816 


1 


IS 


11.5 


14.892 


34.000 


0 . 50 


ioo. a 


-.02546 


- 29.699 


1 


16 


11 . 5 


14. 875 


34 . 000 


0 . 50 


100. 8 


-.03118 


- 36.504 


1 


17 


11.5 


14.858 


34.000 


0.50 


100.8 


- .03477 


- 40. 743 


1 


18 


11.5 


14.837 


34.000 


0.50 


100.8 


-.03600 


- 42.198 


1 


19 


11.6 


14.820 


34.000 


0.50 


100.8 


-.03477 


- 40.904 


1 


20 


11.7 


14.304 


34.000 


0.50 


100.8 


-.03118 


- 36.911 


1 


21 


11.6 


14.787 


34.000 


0.50 


100.8 


-.02546 


- 29.966 


1 


22 


11 . 7 


14. 770 


34.000 


0.50 


100.8 


-.01800 


- 21.230 


1 


23 


11.9 


14.754 


34.000 


0.50 


100 . 8 


-.00932 


- 11.108 


1 


24 


12.1 


14.737 


34.000 


0.50 


100.8 


0.0 


0.0 




1 


12 . 4 


14. 728 


34.000 


0.50 


- 3.5 


0.0 


0.0 


2 


2 


12.5 


14 .727 


34.000 


0.50 


- 100.8 


0.0 


0.0 


2 


3 


12.7 


14.731 


34.000 


0.50 


- 134. 3 


0.0 


0.0 


2 


4 


12.3 


14.741 


34.000 


0.50 


- 248.4 


0.0 


0.0 


2 


5 


12.9 


14. 753 


34.000 


0.50 


- 288.7 


0.0 


0.0 


2 


6 


13.0 


14.765 


34.000 


0.50 


- 302.5 


0.0 


0.0 


2 


7 


13.1 


14.777 


34.000 


0.50 


- 288.9 


0.0 


0.0 


2 


8 


13.2 


14. 786 


34.000 


0.50 


- 248. 7 


0.0 


0.0 


2 


9 


13.3 


14.791 


34.000 


0.50 


- 184.7 


0.0 


0.0 


2 


10 


13.4 


14.790 


34.000 


0.50 


- 101.3 


0.0 


0.0 


2 


11 


13.6 


14.784 


34 . 000 


0.50 


- 4.1 


0.0 


0.0 


2 


12 


13.3 


14.769 


34.000 


0.50 


100.2 


0.0 


0.0 


2 


13 


13 . 9 


14 . 755 


34.000 


0.50 


100.8 


-.00932 


- 12.973 


2 


14 


13.8 


14.742 


34.000 


0.50 


100.8 


-.01800 


- 25.062 


2 


15 


13.6 


14. 729 


34.000 


O.SO 


100.8 


-.02546 


- 35.120 


2 


16 


13.4 


14.716 


34.000 


0.50 


100.8 


-.03118 


- 42. 369 


2 


17 


13.1 


14. 702 


34.000 


0.50 


100.8 


-.03477 


- 46.369 


2 


18 


12 . a 


14.636 


34.000 


0.50 


100 . a 


-.03600 


- 47.034 


2 


19 


12.6 


14.671 


34.000 


0.50 


100 . a 


-.03477 


- 44.547 


2 


20 


12.4 


14.656 


34.000 


0.50 


100.8 


-.03118 


- 39.233 


2 


21 


12.3 


14.641 


34.000 


O.SO 


100.8 


-.02546 


- 31.643 


2 


22 


12.3 


14.625 


34.000 


0.50 


100.8 


-.01800 


- 22.265 


2 


2 3 


12.4 


14.610 


34.000 


0.50 


100.8 


-.00932 


- 11.568 


2 


24 


12.6 


14.590 


34.000 


0.50 


100.8 


0.0 


0.0 
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OAT 


HR 


H 


T 


s 


TAO 


QNET 


WOO 


OH 


1 


1 


49.5 


15.001 


34.000 


0.50 


- 55.7 


0.0 


0.0 


1 


2 


18 . 1 


15.009 


34.000 


0.50 


- 201.5 


0.0 


0.0 


1 


3 


16.2 


15.024 


34.000 


0.50 


- 326.8 


0.0 


0.0 


1 


4 


13.0 


15.047 


34.000 


0.50 


- 423.0 


0.0 


0.0 


1 


5 


12.1 


15.075 


34.000 


0.50 


- 483.5 


0.0 


0.0 


1 


6 


11.8 


15.105 


34.000 


0.50 


- 504.2 


0.0 


0.0 


1 


7 


12.1 


15.132 


34.000 


0.50 


- 483.7 


0.0 


0.0 


1 


a 


12.5 


15.154 


34.000 


0.50 


- 423.5 


0.0 


0.0 


1 


9 


13.0 


15.168 


34.000 


0 . 50 


- 327.5 


0.0 


0.0 


1 


10 


13.5 


IS . 175 


34.000 


0.50 


- 202.4 


0.0 


0.0 


1 


11 


14.2 


15.172 


34.000 


0 . 50 


- 56.6 


0.0 


0.0 


l 


12 


15.0 


15.160 


34.000 


0.50 


99.9 


0.0 


0.0 


1 


13 


15.7 


15.149 


34.000 


0.50 


100.8 


-.00932 


- 14.731 


1 


14 


16.3 


15.138 


34.000 


0 . 50 


100.8 


-.01800 


- 29.577 


1 


15 


16.7 


15.129 


34.000 


0.50 


100.8 


-.02546 


- 43.104 


1 


16 


17.0 


15.119 


34.000 


0.50 


100.8 


-.03118 


- 53.910 


1 


17 


17.3 


15.110 


34 . OOO 


O.SO 


100.8 


-.03477 


- 61.301 


1 


18 


17.6 


15.101 


34.000 


O.SO 


100.8 


-.03600 


- 64.690 


1 


19 


13.0 


15.092 


34.000 


0.50 


100.8 


-.03477 


- 63.643 


1 


20 


18.4 


15.084 


34.000 


0.50 


100.8 


-.03118 


- 58.407 


1 


21 


19.1 


15.076 


34.000 


0.50 


100.8 


-.02546 


- 49.130 


1 


22 


19.9 


15.068 


34.000 


0.50 


100.8 


-.01800 


- 36 . 12 S 


1 


23 


20.9 


15.061 


34.000 


O.SO 


100.8 


-.00932 


- 19.604 


1 


24 


22 . 3 


15.054 


34.000 


0.50 


100.8 


-.00000 


- 00.001 


2 


1 


23.1 


15.054 


34.000 


O.SO 


- 55.7 


0.0 


0.0 


2 


2 


23.1 


15.061 


34.000 


0.50 


- 201. 5 


0.0 


0.0 


2 


3 


14.9 


15.077 


34.000 


3.50 


- 326.8 


0.0 


0.0 


2 


4 


13.3 


15.100 


34.000 


0.50 


- 423.0 


0.0 


0.0 


2 


5 


12.0 


15.128 


34.000 


0.50 


- 483.5 


0.0 


0.0 


2 


6 


11.7 


15.157 


34.000 


0.50 


- 504.2 


0.0 


0.0 


2 


7 


12 . 1 


15.134 


34.000 


0 . 50 


- 483.5 


0.0 


0.0 


2 


a 


12.4 


15.206 


34.000 


0.50 


- 423.5 


0.0 


0.0 


2 


9 


12.9 


15.221 


34.000 


0.50 


- 327.5 


0.0 


0.0 


2 


10 


13.5 


15.227 


34.000 


0.50 


- 202.4 


0.0 


0.0 


2 


11 


14.2 


15.224 


34.000 


0.50 


- 5 6.6 


0.0 


0.0 


2 


12 


15.0 


15.212 


34.000 


0.50 


99.9 


0.0 


0.0 


2 


13 


15.6 


15.201 


34.000 


O.SO 


100.8 


-.00932 


- 14.641 


2 


14 


16.2 


15 . 191 


34.000 


0.50 


100.8 


-.01800 


- 29 . 359 


2 


IS 


16.6 


15.181 


34.000 


0.50 


100.8 


-.02546 


- 42 . 771 



16 . ) 
17.2 

17.6 
17 . > 

ia .2 

18.7 
19.1 

19.7 
20.4 



15.172 

15.162 

15.153 

15.144 

15.136 

15.128 

15.119 

15.112 

15.104 



34.000 

34.000 

34.000 

34.000 

34.000 

34.000 

34.000 

34.000 

34.000 



0.50 

0.50 

0.50 

0.50 

0.50 

0.50 

0.50 

0.50 

0.50 



100.8 

100.8 

100.8 

100.8 

100.8 

100.8 

100.8 

100.8 

100.8 



-.03118 

-.03477 

-.03600 

-.03477 

-.03118 

-.02546 

-.01800 

-.00932 

-.00000 



- 53.606 
- 61 .038 
- 64 .530 
- 63.417 
- 57.717 
- 48.092 
- 34.610 
- 18. 398 
- 00.001 



(C) 



Figure :6. (a) and (b) same as figure (1). except for 

8 -- -- ’<» (c) and (d) same as (a) 



Case lie with peak at hour 18. 
and (b), except for Case lid. 
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QHET 
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OH 


1 


1 


6.3 


14. 9S6 


34.000 


3.50 


- 3.5 


0.0 


0.0 


1 


2 


7.2 


14.939 


34.000 


0.50 


-100.8 


- .00932 


- 6.7S8 


1 


3 


7.3 


14.935 


34.000 


0. SO 


-184 . 3 


-.01800 


-14.212 


1 


4 


8. 1 


14.939 


34.000 


0.50 


-248.4 


-.02546 


-20.986 


1 


S 


8. 3 


14.947 


34.000 


0. SO 


-288. 7 


-.03118 


-26.310 


1 


6 


8. 4 


14.9S4 


34.000 


0.50 


-302 . S 


-.03477 


-29.781 


1 


7 


8.S 


14.963 


34.000 


0.50 


-288.9 


-.03600 


-31.198 


1 


8 


8 . 6 


14.969 


34.000 


0.50 


-248. 7 


-.03477 


-30.S26 


1 


9 


8.8 


14.969 


34.000 


0.50 


-184. 7 


-.03118 


-27.771 


1 


10 


9.8 


14 . 962 


34.000 


0.50 


-101. 3 


- . 02 S 46 


-22.787 


1 


11 


9.0 


14.946 


34.000 


0.50 


- 4. 1 


-.01800 


-16. 328 


1 


12 


9.2 


14.920 


34.000 


0.50 


*100.2 


-.00932 


- 8.6S2 


1 


13 


9.6 


14.896 


34.000 


0.50 


100.8 


0.0 


0.0 


1 


14 


10.0 


14 .874 


34.000 


3.50 


100.8 


0.00932 


9.319 


l 


IS 


10.6 


14.8S3 


34.000 


0.50 


100.8 


0.01800 


18.919 


1 


16 


11.3 


14.835 


34.000 


0. SO 


100.8 


0.02S46 


28.302 


1 


17 


12.0 


14.817 


34.000 


0.50 


100.8 


0.03118 


36.838 


1 


18 


12.8 


14.797 


34.000 


0.50 


100.8 


0.03477 


43. 781 


1 


19 


13.6 


14 .783 


34.000 


O.SO 


100.8 


0.03600 


48.172 


1 


20 


14.4 


14. 770 


34.000 


0.50 


100.8 


0.03477 


49.208 


1 


21 


15.1 


14.758 


34.000 


O.SO 


100.8 


0.03118 


46.482 


1 


22 


IS . 8 


14.747 


34.000 


O.SO 


100.8 


0.02S46 


39.7SS 


l 


23 


16.4 


14.736 


34.000 


0.50 


100.8 


0.01800 


29. 195 


1 


24 


16. 7 


14. 72S 


34.000 


O.SO 


100.8 


0.00932 


1S.S26 


2 


1 


16.9 


14.720 


34.000 


0.50 


- 3.5 


0.0 


0.0 


2 


2 


16.9 


14.722 


34.000 


O.SO 


-100.3 


-.00932 


-IS. 818 


2 


3 


16-7 


14 . 728 


34. 000 


O.SO 


-184.3 


-.01800 


-30.300 


2 


4 


16. 3 


14 .737 


34.000 


0.50 


-248.4 


-.02546 


-42.078 


2 


5 


IS. 9 


14 . 749 


34.000 


0. 50 


-288.7 


-.03118 


-SO. 22 3 


2 


6 


IS. 4 


14 . 761 


34.000 


0.50 


-302.5 


-.03477 


-54. 338 


2 


7 


14.9 


14.773 


34.000 


O.SO 


-288.9 


-.03600 


-S4 . 494 


2 


3 


14.4 


14. 783 


34.000 


0.50 


-248. 7 


-.03477 


-S1.081 


2 


9 


14.1 


14 . 788 


34.000 


0.50 


-184. 7 


-.03118 


-44.649 


2 


10 


13.9 


14 . 788 


34.000 


O.SO 


-101.3 


-.02546 


-35.798 


2 


11 


13.8 


14.782 


34.000 


O.SO 


- 4.1 


-.01800 


-2S.074 


2 


12 


13.9 


14 . 768 


34 . 000 


0.50 


100.2 


-.00932 


-12.981 


2 


13 


14.1 


14.7SS 


34.000 


0.50 


100.8 


0.0 


0.0 


2 


14 


14. 4 


14.742 


34.000 


0. SO 


100.8 


0.00932 


13.323 


2 


IS 


14.8 


14.730 


34. 000 


O.SO 


100.8 


0.01300 


26.414 


2 


16 


15.4 


14.718 


34.000 


0. 50 


100.8 


0.02546 


38.631 


2 


17 


16.0 


14 . 70 7 


34.000 


O.SO 


100.8 


0.03118 


49.224 


2 


18 


16.8 


14.695 


34.000 


0. 50 


100.8 


0.03477 


S7. 351 


2 


19 


17.6 


14.685 


34.000 


0.50 


100.8 


0.03600 


62.131 


2 


20 


18.4 


14.675 


34.000 


O.SO 


100.8 


0.03477 


62.760 


2 


21 


19.0 


14.667 


34.300 


O.SO 


100.8 


0.03118 


58.667 


2 


22 


19.8 


14.659 


34.000 


0.50 


100.8 


0.02S46 


49.687 


2 


23 


20.3 


14.651 


34.000 


O.SO 


100.8 


0. 01800 


36. 190 


2 


24 


20.6 


14.642 


34.000 


0.50 


100.8 


0.00932 


19.137 



(a) 




1 12 21 

TIME (HR) 

(b) 




12. CO : 2 - SO 13.00 ’.3.50 14.00 14.50 15.00 15.50 16. CO IS. SO 

“"EMP . (RELATIVE TO - I RST) 



TOO 
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OAf 


SR 


a 


T 


S 


TAU 


QNET 


woo 


OH 


1 


X 


6 . 0 


14.956 


34.000 


0.50 


- 3.5 


- .03600 


-22.140 


1 


2 


6. 9 


14.937 


34.000 


0.50 


-100. 8 


- .03477 


-24.417 


1 


3 


7.3 


14.932 


34.000 


0.50 


-184.3 


-.03118 


-23.270 


1 


4 


7.8 


14.934 


34.000 


0.50 


-248.4 


-.02546 


-20.008 


l 


5 


a.x 


14.942 


34.000 


0.50 


-288.7 


-.01800 


-14.652 


: 


6 


3.3 


14.949 


34.000 


0.50 


-302.5 


-.009 32 


- 7.808 


i 


7 


3. 7 


14. 959 


34.000 


O.SO 


-288.9 


0.0 


0.0 


i 


8 


9 . X 


14.965 


34 .000 


0.50 


-248. 7 


0. 009 32 


8.440 


l 


9 


9.5 


14.966 


34.000 


0.50 


-184.7 


0.01800 


17.029 


l 


10 


XO. X 


14.960 


34.000 


0 . 50 


-101.3 


0.02546 


25.370 


i 


XX 


10.7 


14.948 


34.000 


0.50 


- 4.1 


0.03118 


32.978 


X 


12 


XX. 5 


14.928 


34.000 


0.50 


100.2 


0.03477 


39.220 


l 


13 


12.2 


14.91X 


34.000 


O.SO 


100.8 


0.03600 


43.180 


i 


14 


X3.0 


14. 896 


34.000 


0.50 


XOO. 8 


0.03477 


44.276 


i 


IS 


13.7 


14.382 


34.000 


0.50 


100. 8 


0.03118 


41.902 


X 


16 


X4 . 3 


14.869 


34.000 


0.50 


100.8 


0.02546 


35 . 893 


1 


17 


14.3 


X4. 857 


34. 000 


0.50 


100.8 


0.01800 


26.426 


X 


18 


15.2 


14.844 


34. 000 


0.50 


100.8 


0.00932 


14. Ill 


X 


19 


15 .5 


14.833 


34.000 


0.50 


100.8 


0.0 


0.0 


L 


20 


15.6 


14 . 822 


34.000 


0. 50 


100.8 


-.00932 


-14.562 


X 


2 X 


15.5 


14.8X1 


34.000 


0.50 


100.8 


-.01800 


-28.200 


X 


22 


15.4 


14.800 


34.000 


0.50 


xoo. a 


-.02546 


-39.680 


1 


23 


15 . 2 


14 . 789 


34.000 


0.50 


XOO. 8 


-.03118 


-48.077 


X 


24 


14.9 


14.776 


34.000 


O.SO 


100.8 


-.03477 


-52.860 


2 


X 


14.6 


X 4 . 771 


34.000 


O.SO 


- 3.5 


-.03600 


-53.604 


2 


2 


14.3 


14.771 


34.000 


0.50 


-100.9 


-.03477 


-50.584 


2 


3 


14.0 


14 . 776 


34.000 


0.50 


-184.3 


-.03118 


-44.329 


2 


4 


13.7 


14.786 


34.000 


0.50 


-248.4 


-.02546 


-35. 447 


2 


S 


X 3 . 6 


14. 798 


34.000 


0.50 


-288. 7 


-.01800 


-24.685 


2 


6 


13.5 


14.8X0 


34.000 


0.50 


-302.5 


-.00932 


-12.683 


2 


7 


13.6 


14.822 


34.000 


0.50 


-288.9 


0.0 


0.0 


2 


a 


X 3 . 9 


14.831 


34.000 


0.50 


-248. 7 


0.009 32 


12.864 


2 


9 


14.2 


14.837 


34.000 


0 . 50 


- X 84 . 7 


0.01800 


25.415 


2 


xo 


X 4 . 8 


X4. 837 


34.000 


0.50 


-XOX. 3 


0.02546 


37.112 


2 


XX 


XS. 4 


14.831 


34.000 


0.50 


- 4.X 


0.03118 


47.309 


2 


12 


16.2 


14. 3X8 


34.000 


0.50 


100.2 


0.03477 


55.256 


2 


13 


17.0 


14.808 


34.000 


0.50 


100.8 


0.03600 


60 . 00 2 


2 


14 


17.3 


14.798 


34.000 


O.SO 


LOO. 3 


0.03477 


60.692 


2 


IS 


18. 5 


14.799 


34.000 


0.50 


100.8 


0.03118 


56.796 


2 


16 


19 . 2 


14.781 


34. 000 


0.50 


100.8 


0.02546 


48.153 


2 


17 


19.7 


X 4 . 7 7 3 


34.000 


0. 50 


100.3 


0.01800 


35.109 


2 


IS 


20.0 


14 . 764 


34.000 


0.50 


100.8 


0.00932 


18.583 


2 


19 


20.2 


14. 756 


34 .000 


0.50 


100.3 


0.0 


0.0 


2 


20 


20.2 


14 . 749 


34.000 


0.50 


loo . a 


- . 00932 


-18.903 


2 


2 X 


20.0 


14. 741 


34.000 


0.50 


100.3 


-.01800 


-36.339 


2 


22 


19 . 7 


X4. 733 


34.000 


0. 50 


loo . a 


-.02546 


-50.746 


2 


23 


19.3 


X4. 725 


34.000 


0.50 


100.8 


-.03118 


-60.996 


2 


24 


18. 8 


14 .7X5 


34.000 


0.50 


100.3 


-.03477 


-66.499 
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10.8 
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IS. 7 

16.2 

16.7 
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■3 14.9430 
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14.731 
14. 720 
14.709 
14.699 
14.689 
14.679 
14.670 



34.000 
34.000 
34.000 
34.000 
34.000 
34.000 
34.000 
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34.000 
34.000 
34.000 
34. 000 
34.000 
34.000 
34.000 
34.000 
34.000 
34.000 
34.000 
34.000 
34.000 
34.000 
34.000 
34.000 

34.000 
34.000 
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Figure 30. Same as figure (28), 
for reversel of peak directions. 
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Figure 31. Contours as in figure (1), except for Case Ilia in two dimensional mode, 
fa), (b) , fc) , and (d) are hour 6, 12, 18 and 24 respectively. 
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Figure 37. Same as figure (3S), except for 
reversel of peaks. 
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Figure 38. Sane as figure (31), except for Case Illb. 
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Figure 41. Contours as in figure (1), except for Case IVa. (a), (b) , (c) , and (d) have the maximum 
upward notion at hour 6, 24, 18 and 12 respectively. 
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Figure 44. Same as figure (31), except for Case IVa. 
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Figure 46. Same is figure (33), 
except for Case IVa. 



119 



"1 



o 

3 

co 



a 

o 

J 7 

X 



15.0000 
15.0090 
ts.caio 
:s.ox7o 
15.0750 
15.1050 
15. 1330 
IS. 1550 
15.1690 
15.1750 
is. ; 7*o 
15.1620 
IS. 1500 
15.1300 
15. 12’0 
15.1170 
15. 1090 
15.1000 
15. 09*0 
15. 3960 
15.0930 
15.0790 
15.0730 
15.0670 




-Ad 



13.00 



o 




' 4 . CO 



14.50 15.00 15.50 16.00 

T £* p . (RELATIVE '0 FIRST) 





120 



1 -r-i- ■!— ■ 1 1 i ' > ■ ■■ ; i i i i r i ■ ^ 




Figure 48. Same as figure (31), except for Case IVb. 
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Figure 51. Same as figure (31), except for Case Va. 
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